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EVALUATION  STATEMENT 


This  report  describes  a series  o1  studies  performed  by  Syracuse 
' University  for  Air  Force  Cambridge  Research  Laboratories  and  then  Rome 

Air  Development  Center  under  Contract  No.  FI9628-76-C-0300.  Two  sections 
makf  up  tiie  report.  The  final  section  is  a review  of  the  earlier  scienti- 
fic  'eports  on  the  application  of  the  general  formulation  of  matrix  solu- 
[ tion  to  multiple  regions  and  multiple  sources.  These  general  solutions 

I were  applied  to  the  specific  problems  of  arrays  of  cavity-backed  and  waveguide- 

fed  apertures.  The  first  section  is  a scientific  report  on  the  problem  of 
I plane  wive  scattering  by  a body  of  revolution  of  a homogeneous  mattrial, 

r 

I The  use  of  both  E field  and  H field  combined  and  formulated  in  terms  of  the 

; equivalent  electric  and  magnetic  currents  over  the  body  surface  gives  unique 

I solutions  even  at  body  resonant  frequencies.  The  technique  is  applied  to 

examples  (dielectric  spheres)  and  the  results  compared  to  thos  of  Mfller 
who  uses  a similar  but  different  approach.  Computer  programs  are  a so 


listed.  These  programs  may  have  application  to  Air  Force  programs  uch  as 
missile  plume  identification  and  human  intruder  detection. 
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PART  ONE 


ELECTROMAGNETIC  SCATTERING  FROM  A HOMOGENEOUS  MATERIAL 
BODY  OF  REVOLUTION 

THEORY  AND  EXAMPLES 


I . INTRODUCTION 

Th<?  problem  of  plane-wave  scattering  by  a homogeneous  material  body 
of  revolution  is  formulated  in  terms  of  equivalent  electric  and  magnetic 
currents  over  the  body  surface.  Application  of  boundary  conditions  leads 
to  a set  of  four  integral  equation?  to  be  satisfied.  Linear  ccmbinatlons 
of  these  four  equations  lead  to  a oupled  pair  of  equations  to  be  solved. 
One  choice  of  combination  constants  gives  the  formulation  described  by 
Poggio  and  Miller  [1].  This  formulation  has  been  applied  to  material 
cylinders  by  Chang  and  Harrington  [2],  and  to  material  bodies  of  revolu- 
tion by  Wu  [3].  We  will  call  this  choice  the  PMCHW  formulation  (formed 
by  the  initials  of  the  above  cited  investigators) . 

Another  choici  of  combination  constants  gives  the  formulation 
obtained  by  Muller  [<♦].  This  formulation  has  been  applied  to  dieleci  ric 
cylinders  by  Solodukhov  and  Vasil'ev  [5]  and  by  Morlta  [6],  and  to  bodies 


[1]  A.  J.  Poggio  and  E.  K.  Miller,  "Integral  Equation  Solutions  of  Three- 
dimensional  Scattering  Problems,"  Chap.  4 of  Computer  Techniques  for 
Electromagnetics,  edited  by  R.  Mittra,  Pergamon  Press,  1973,  Equa- 
tion (4.17). 

[2]  Yu  Chang  and  1.  F.  Harrington,  "A  Surface  Formulation  for  Characteristl 
Modes  of  Material  Bodies,"  Report  TR-74-7,  Dept,  of  Electrical  and  Com- 
puter Engineering,  Syracuse  University,  Syracuse,  N.Y.,  October  1974. 

[3]  T.  K.  Wu,  "Electromagnetic  Scattering  from  Arbitrarily-Shaped  Lossy 
Dielectric  Bodies,"  Ph.D.  Dlssertarion,  University  of  Mississippi,  1976. 

[4]  C.  Muller,  Foundations  of  the  Mathematical  Theory  of  Electromagnetic 
Waves,  Springer-Verlag,  1969,  p.  301,  Equations  (40)-(41).  (There  are 
some  sign  errors  in  these  equations.) 

[5]  V.  V.  Solodukhov  and  E.  N.  Vasil'ev,  "Diffraction  of  a Plane  Electromagnetic 
Wave  by  a Dielectric  Cylinder  of  Arbitrary  Cross  Section,"  Soviet  Physics  - 
Technical  Physics,  vol.  15,  No.  1,  July  1970,  pp.  32-36. 

[6]  N.  Morlta,  "Analysis  of  Scattering  by  a Dielectric  Rectangular  Cylinder  by 
Means  of  Integral  Equation  Formulation,"  Electronics  and  Comnunl cat ions  in 
Japan,  vol.  57-B,  No.  10,  October  1974,  pp.  72-80. 


1 


k 


of  revolut  on  by  Vasil'ev  and  Materikova  [7].  We  will  call  this  choice 
the  Mllller  formulation.  Conditions  for  the  uniqueness  of  solutions  are 
established  in  tenrn  of  the  combination  constants.  It  is  found  thai 
solutions  to  both  the  PMCHW  formulation  and  to  Muller's  formulation  are 
unique  at  <11  frequencies. 

Numerical  solutions  to  the  coupled  pair  if  equations  are  obtained  by 
the  method  ( f moments  [8].  It  is  relatively  easy  to  obtain  numerical  sc  Lu- 
tions  to  these  equations  because  the  required  operators  are  the  same  as  chose 
evaluated  in  earlier  reports  [9,  10].  An  exemplary  computer  program  ca  able 
of  obtaining  the  solution  to  both  the  PMCHW  formulation  and  the  TlUller  ormu- 
lation  is  described  and  listed.  This  is  a main  program  which  uses  subr  lutines 
similar  to  those  in  [10]  to  compute  the  equivalent  electric  and  nagnetic 
currents  and  the  two  principal  plane  scattering  patterns  for  a Toss-free 
homogeneous  body  of  revolution  excited  by  an  axially  incident  electromagnetic 
plane  wave.  Computed  results  for  the  equivalent  currents  and  principal  plane 
scattering  patterns  of  a dielectric  sphere  whose  relative  dielectric  const mt 
is  four  show  reasonable  agreement  between  our  solution  to  the  PMCHW  formula- 
tion, our  solutloi  to  the  Muller  formulation,  and  the  "exact"  series  [11] 
solution  in  the  resonance  region.  Computer  program  subroutines  which  calcu- 
late tl  e "exact"  series  solution  for  perfectly  conducting  spheres  as  well  as 
for  loss-free  homogeneous  spheres  will  be  described  and  listed  in  a subsequent 
report . 

[7]  E.  N.  Vasil'ev  and  L.  B.  Materikova,  "Excitation  of  Dielectric  Bodies  of 
Revolution,"  Soviet  Physics  - Technical  Physics,  vol.  10,  No.  10,  April 
1966,  pp.  I401-1A06. 

[8]  R.  F Harrington,  Field  Computation  by  Moment  Metho  .s,  Macmillan  Co., 

New  York,  1968. 

[9]  J.  R.  Mautz  and  R.  F.  Harrington,  "H-Fleld,  E-Fleld,  and  (:ombined  Field 
Solutions  for  Bodies  of  Revolution,"  Interim  Technical  Renort  RADC-TR- 
77-109,  Rome  Air  Development  Center,  Grlfflss  Air  Force  lase.  New  York, 
March  1977. 

[10]  J.  R.  Mautz  and  R.  F.  Harrington,  "Computer  Programs  for  H-Fie^d, 

E-Fleld,  and  Combined  Field  Solutions  for  Bodies  of  Revolution." 

Interim  Technical  Report  RADC-TR-77-215,  Rome  Air  Development 
Center,  Grifflss  Air  Force  Base,  New  York,  June  1977. 

[11]  R.  F.  Hairington,  Time- Harmonic  Electromagnetic  Fields,  McGraw-Hill 
Book,  Co.,  1961.  Section  6-9. 
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II.  SURFACE  INTEGRAL  EQUATION  FORMULATION 

An  electromagnetic  field  propagating  In  a homogeneous  medluM  of 
permeability  and  permittivity  is  incident  on  the  surface  Sofa 
homogeneous  obstacle  of  permeability  and  permittivity  The  sub- 

script e denotes  exterior  medium  and  the  subscript  d denotes  diffractln  ; 
medium.  We  wish  to  calculate  the  scattered  electromagnetic  field  H ' 
outside  S and  the  diffracted  electromagnetic  field  E,  H inside  in  terms 
of  the  electromagnetic  field  E^^,  which  would  exist  on  S in  t le  absence 
of  the  obstacle.  This  original  problem  is  shown  in  Fig.  1 whera  J^, 
are  the  electric  and  magnetic  sources  of  E^,  and  n is  the  ui  it  normal 
vector  which  points  outward  from  S. 


The  equivalence  principle  (stated  in  Append* x A)  is  used  to  piece 

s s 

together  an  outside  situation  cons  Lsting  of  medium  p , e and  field  E , H 

0 0 ^ .AA' 

outs  Lde  S and  an  inside  situation  consisting  of  me  Hum  p , e and  field 
1 i e e 

-E  , -H^  li  side  S.  Tl.is  composite  situation  is  shown  in  Fig.  2.  Since 

s s i i 

E , II  is  source-free  outside  S and  E , H is  source-free  inside  S,  the 

Only  sources  in  Fig.  2 are  the  equivalent  electric  surface  current  J and 

the  equivalent  magnetic  surface  current  M on  S. 


As  a second  application  of  the  equivalence  principle,  we  combine  an 
outside  situation  consisting  of  medium  p^,  and  zero  field  with  an  li  side 
situation  consisting  of  medium  p^,  i ^ and  field  E,  H.  This  combinatio  i of 
situations  is  shown  in  Fig.  3.  Sln< e E,  H is  source-free  inside  S,  tie  only 
sources  n Fig.  3 are  the  equivalent  electric  surface  current  -J  and  the 
equivalent  magnetic  surface  current  -M  on  S.  By  using  (A-1)  and  (A-2)  to 
express  the  surface  currents  in  terms  of  the  discontinuities  of  the  tangen- 
tial fields  across  S and  by  using 


n X E = n X (E®  + E^) 


n X !l  = n > (H®  + H^) 


(1) 

(2) 


on  S,  the  Interested  reader  can  verify  that  the  surfs'ie  currents  in  Fig.  3 
are  Indeed  the  negatives  of  those  in  Fig.  2.  Equations  (1)  and  (2)  are  the 
boundary  conditions  that  the  tangential  components  of  the  field  i in  the 
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original  problem  as  shown  in  Fig.  1 are  continuous  across  S. 


I 


The  scattered  field  E®,  outside  S and  the  diffracted  f eld  E,  H_ 
inside  S could  easily  be  calculated  if  J and  M were  known  becaus  ; the  media 
into  which  J and  M radiate  is  homogeneous  in  Figs.  2 and  3.  We  nave  to 
determine  J and  M.  The  equivalence  principle  rtates  that  there  exist  J and 
M which  radiate  the  fields  in  Figs.  2 and  3,  but  the  equivalence,  principle 
does  not  tell  what  J and  M are.  The  equlva  ence  principle  does  state  that 


c> 

It 

>->  i 

X H 

(3) 

M = E 

X n 

(4) 

A*.  A.*- 

but  this  is 

not  very  useful 

because  E and  H are  unknown. 

From 

Figs.  2 and  3, 

i 

(>) 

-n 

X E = 

n X E 

t*-e 

— 

i 

o) 

-n 

X u"  = 

n X n 

-n 

II 

+ y 

X 

0 

(7) 

-n 

II 

X 

0 

(8) 

where 


E is  the  electric  field  just  inside  S due  to  J,  M,  radiating  in  u ,e 
0*^  ^ ^ 6 6 

H is  the  magnetic  field  ju5t  inside  S due  to  J,  M,  radiating  in  p ,c 

e e 

Ej  is  the  elei  trie  field  juit  outside  S due  to  J,  M,  radiating  in  P,,g. 

•~a  da 

is  the  magn  !tlc  field  juit  outside  S due  to  J,  M,  radiating  rn 

The  equivalent  currents  J,  M which  appear  in  Figs.  2 and  3 satisfy  (5) -(8) 
because  (5) -(8)  were  obtained  from  Figs.  2 and  3.  It  is  shown  in  Appendix  B 
that  the  solution  to  (5) -(8)  is  unique.  Therefore,  (5) -(8)  tmiquely  deter- 
mine the  equivalent  currents  J,  ^ of  Figs.  2 and  3. 

Equations  (5) -(8)  form  a set  of  four  equations  in  the  two  unknowns 
J and  M^.  The  usual  methods  of  equation  solving  apply  only  wht  n the  number 
of  equations  is  equal  to  the  number  of  unknowns.  We  want  to  reduce  the 
set  of  four  equations  (5) -(8)  to  two  c^quatlons.  One  way  to  do  this  is  to 
form  the  1 Lnear  combination 
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(9) 


-r  X (E  + aE^)  = a x 
■“  ->*€  **<1 


of  (5)  and  (7)  and  the  linear  conih  Ination 


-n  X (H  + BH't)  = n < 

of  (6)  and  (8)  where  a and  6 are  complex  constants. 

The  solution  J,  M to  (5)-(8)  satisfies  (9)  and  (10).  This  J,  M 
will  be  the  only  solution  to  the  pair  of  equations  (9)  and  (10)  if 


-n  X (E  + aE";)  = 0 


-n  X (H  + t;H^)  = 0 
*-€  ■'.-d 

have  only  tha  trivial  solution  J = M = 0.  From  (11)  and  (12), 


P = - a6  P,  (13) 

e a 

where  P is  the  complex  power  flow  of  E , H inside  S ani  P,  is  the  complex 
e d 

power  flow  of  Id’  5d  outside  S.  Tie  asterisk  in  (13)  denotes  complex  con- 
jugate. If  a6  is  real,  then  the  real  part  of  (13)  reduces  to 

Real(P^)  = - a6*Real(P,)  (14) 

e a 

* 

If  a0  is  not  only  real  but  also  positive,  then 


Rea:(P_,)  = 0 
d 

because  both  Real(P^)  and  Rc*al(P^)  are  greater  than  or  equal  to  zero. 


Equation  (15)  implies  that 


nxEj=nxH^=0 


Substitution  of  (16)  into  (11)  and  (12)  yields 


nxE  = nxH  = 0 
— — e — -e 
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The  system  of  equations  (16)  and  (17)  Is  preclsel  / the  homogeneous  system 

of  equations  associated  with  (5)-(8).  It  w.as  shovn  In  Appendix  B that 

this  homogeneo  IS  system  of  equations  las  only  the  trivial  solutloi  J = M = 0 
* ^ 

Therefore,  If  t6  Is  real  and  posltlv  then  the  coupled  pair  of  < quatlons 

(11)  and  (12)  has  only  the  trivial  solution  J = M = 0 so  that  the  solution 

M to  (5)-(8)  Is  the  only  solution  to  the  coupled  pair  of  equations  (9) 

and  (10). 

If  a = 8 = 1,  then  (9)  a id  (10)  become 

-n  X (e“  + Ej)  = n X (18) 

-n  X (h“  + n't)  = n X 19) 

The  set  of  equati ms  (18)  and  (19)  Is  the  coupled  pair  of  surface  Integral  | 

equations  transcr  bed  by  Pogglo  and  Miller  [1].  We  call  these  equations  | 

the  PMCHW  equations.  Since  a = 6 = 1 Implies  that  xB  Is  real  and  positive,  ] 

the  argument  cons  sting  of  (11) -(17)  and  Involving  real  power  flow  shows  i 

that  (18)  and  (19)  uniquely  determine  the  desired  J,  M of  Figs.  2 and  3.  • 

That  (18)  and  (19)  uniquely  determine  M of  Figs.  2 and  3 can 
also  be  shown  as  follows.  The  desire  1 J,  M of  Figs.  2 and  3 satisfies 

(18)  and  (19)  because  (18)  and  (19)  were  obtained  from  Figs.  2 and  3. 

This  desired  J,  M will  be  the  only  solution  to  (18)  and  (19)  if  the  assocl- 
ated  set  of  liomogeneous  equations 


-n  X (e“  + Ej)  = 0 

(20) 

-n  X (h  + Hj)  = 0 
— — e —d 

(21) 

has  only  the  trivial  solution  J = M - 0. 

/ The  following  argument  shows  that  (20)  and  (21)  have  only  the  trivial 

solution  d " M = 0.  Let  E^,  be  the  electromagnetic  field  outside  S due 

to  J,  M radiating  in  e^.  Let  be  the  electromagnetic  field  Inside 

S due  to  J,  H radiating  in  e^.  Use  the  equivalence  principle  to  form 

the  composite  situation  consisting  of  medium  p^,  and  field  E j,  outsit  3 

S and  medium  p , e„  and  field  -E  , -H  inside  S as  shown  In  Fig.  4.  In 

Fig.  A,  Ej,  Hj  Is  a source-free  Maxwellian  field  outside  S.  Since  E , H is 
**-Q  •*'0  —e 
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Fig.  4.  Composite  situatior  used  to  prove  that  (20) 
and  (21)  have  only  the  trivial  solution 
J = M = 0. 


1 ig.  5.  Composite  situation  used  to  prove  that  (26) 
and  (27)  have  only  the  trivial  solution 
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a source-free  Maxwellian  field  Inside  S,  the  field  appearing  In 

Fig.  4 Is  also  a source-free  Maxwellian  field  Inside  S.  Now,  (20)  and  (21) 
state  that  the  tangential  components  of  the  field  in  Fig.  4 are  continuous 
across  S.  Thus,  Fig.  4 Is  entirely  source-free  so  that  the  field  In  Fig.  4 
is  zero  everywhere  In  which  case  (B-L)-(B-4)  are  satisfied.  But,  as  shown  i.i 
Appendix  B,  (B-l)-(B-4)  have  only  the  trivial  solution  J = M = 0.  Hence, 

(20)  and  (21)  have  only  the  trivial  solution  J = M « 0. 

If 

0.  = - ^ 


(22) 

^23) 


then  (9)  and  (10)  become 


-n  X (E ^ E_,)  = n X E 

- .-e  e -d 

e 


(24) 


(H  - H .) 
-.e  p --d 
e 


n X H 


(25) 


The  set  of  e4uatlons  (24)  and  (25)  Is  the  coipled  pair  of  surface  integral 

equations  obtained  by  Muller  [4].  We  ;all  t lese  equations  the  Muller  eq  la- 

tlons.  The  singularity  that  the  kernels  of  he  integral  equations  (24) 

and  (25)  exhibit  as  the  source  point  passes  through  the  field  point  is  lot 

as  pronounced  as  the  singularity  of  the  kemsls  of  (18)  and  (]9).  If 

* 

Eg,  pj  , and  e^  are  real  In  (22) -(25),  then  “6  is  real  and  positive. 

In  this  case,  the  argument  consisting  of  (11) -(17)  shows  that  (24)  and  (25) 
uniquely  determine  the  desired  J,  M of  Figs.  2 and  3. 

An  alternate  proof,  valid  for  lossy  media,  that  (24)  ind  (25)  uniqu  ly 
determine  the  desired  J,  M is  presented.  This  proof  is  similar  to  the  argu- 
ment which  used  Fig.  4 to  show  that  (18)  and  (19)  uniquely  determine  the 
desired  J,  M and  is  as  follows.  The  desired  J,  M of  Figs.  2 and  3 satisfies 
(24)  and  (25)  because  (24)  and  (25)  were  obtained  from  Figs.  2 and  3.  This 
desired^,  ^ will  be  the  only  solution  to  (24)  and  (25)  if  the  associ;  ';d 
set  of  homogeneous  equations 
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(26) 


-n  x (e"  - — Et) 
e 


-n  X (jT  - 7^  H^)  - 0 
— ->-€  y 

e 


haa  only  the  trivial  solution  J = M 


The  following  argument  shows  that  (26)  and  (27)  have  only  the  trivial 

solution  J *=  M = 0.  Let  be  the  electromagnetic  field  outside  S due 

to  _J,  ]M  radiating  in  u^,  e^.  Because  the  electromagnetic  field  is 

a source-free  Maxwellian  field  in  outside  S,  the  dual  electromagnetic 

field  n . H ,»  — ^ E , where 
d rij 

d 


n 


Is  also  a source-free  Maxwellian  field  outside  S.  Let  E , H be  the 

-e  >“€ 

electromagnetic  field  inside  S due  to  J,  M radiating  in  p , e . Because 

— w>.  e e 

the  electromagnetic  field  E , H is  a source-free  Maxwe]llan  field  inside  S, 

--e  — e . 

the  dual  electromagnetic  field  n H , E where 

e n 

e 


is  also  a source-free  Maxwellian  field  inside  S.  Use  the  equivalence 

principle  to  form  the  composite  situation  consisting  of  medium  p,,  e.  and 
^ d d 

field  n . H » ^ outside  S and  medium  p , e and  field 


u e 

(n  H , 


— ^)  Inside  S as  shown  in  Fig.  5.  Now,  (26)  and  (27) 
e 


state  that  the  tangential  components  of  the  field  in  Fig.  5 are  continuous 
acrost  S.  Thus,  Fig.  5 is  entirely  source-free  so  that  the  field  in  Fig.  > 
is  zero  everywhere  in  which  case  (B-l)-(B-4)  are  satisfied.  But,  as  shown 
in  Appendix  B,  (B-l)-(B-4)  have  only  the  trivial  solution  J - M - J.  Hence, 
(26)  and  (27)  have  only  the  trivial  solution  J ■ M - 0. 
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III.  METHOD  OF  MOMENTS  SOLUTION  FOR  BODY  OF  REVOLUTION 


In  this  section,  a method  of  moments  solution  to  (9)  and  (10)  is 
developed  for  a homogeneous  loss-free  body  of  revolution.  Special  cases  of 
(9)  and  (10)  are  the  PMCHW  equations  (18)  and  (19)  and  the  Muller  equations 
(24)  and  (25). 


For  compatibility  with  equation  (40)  on  page  14  of  [9],  we  rewrite 

(9)  as 


+ aEj) 

Mrd  tan 


1 


E^ 

-<an 


(30) 


where  tan  denotes  tangential  components  on  S and  is  given  by  (29).  The 
fields  on  the  left-hand  sides  of  (30)  and  (10)  are  written  as  the  sum  of 
fields  due  to  J and  fields  due  to  M.  Advantage  is  taken  cf  the  fact  that 
the  operator  which  gives  the  electric  i ield  due  to  a magnetic  current  is 
the  negative  of  tha  operator  which  gives  the  magnetic  field  due  to  an 
electric  current  ajid  that  the  operator  which  gives  the  magnetic  field  <lue 
to  a magnetic  current  is  the  square  of  the  reciprocal  of  the  intrinsic 
impedance  times  the  operator  which  gives  the  electric  field  due  to  an 
electric  current.  In  view  of  the  above  considerations,  (30)  and  (10)  be- 
come 


(-  — E (J)  + — H (M)  - — E,(J)  + — Hj(M))^  = — e' 

q~e—  Ti«Ki~  q-d*-  tan  q **  an 
e e e e e 


(31) 


-n  X (h2(J)  + ^ E (M)  + 6Hj(J)  + Ay  E,(M))  = n x 


e - 


(32) 


where  ^ denotes  the  operator  which  gives  the  electric  field  due  to  an 
electric  current.  The  subscript  e or  d on  E denotes  radiation  in  either 
tg  Of  Uj.  e^.  The  superscript  + or  -,  if  present  on  E,  denotes  field 
evaluation  either  just  outside  S or  just  inside  S.  The  H's  in  (31)  and 
(32)  are  the  corresponding  magnetic  field  due  to  electric  current  opeiators. 
We  stress  that  all  E's  and  ^'s  in  (31)  and  (32)  aie,  by  definition,  o erators 
which  give  electric  and  magnetic  fields  due  to  electric  currents,  eve  though 
these  operators  act  on  both  electric  and  magnetic  currents  J and  M in  (31) 
and  (32). 
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Let 


j = I ! d". 

- n=-«  j=l 

(33) 

M = n I I 'v'^.  + V**,  j'*’,) 

"■  e ^ nj  -nj  nj  —nj 

n=-oo  j=l  -j  j j j 

(34) 

where  1^.,  I*^.,  V',,  and  V*^.  are  coe  flclents  to  be  determine  1 and 
nj  nj  nj  nj 

(35) 

1 

f- 

( 

4 ■ 

(16) 

In  (35)  and  (36),  t is  the  arc  length  along  the  generating  curve  of  the  body 

body  of  revolution  and  (J>  isi  the  longitudinal  angle,  u^  and  u,  are  unit 

«-t  **9 

vectors  in  the  t and  ip  directions  respectively  such  that  u.  u^  = n ai  J 
fj(t)  is  the  scalar  function  of  t defined  on  page  10  of  [9].  The  body  of 
revolution  and  coordinate  system  are  shown  in  Fig.  6.  Substitution 
(33)  and  (3A)  into  (31)  and  (32)  yields 


oo 


I 


n=-“ 


tan 


(H'(J\) 

-e  -Mtij 


aHt(j’*’.)),.  V'*’ 

"d  -^nj  tan  nj 


+ 


(_  _e2£nj d . jt  ^ (_:-e  rnl d -d  >nj...  j. 

He  tan  nj  n ^ tan 


^.1 


(37) 


-n  X I 
n»- 


y \ + —B.  rArpj-)  -iiL-)v'<>  + 

j-l  V ’’e  ^d  ^d  % ^d  ^d 


■*•  bh;|;(^,))i^,  + (H  (jj;,) 


+ 6Ht(j‘^ .))!’’’.  I - n X 


'I^e'-nj'  nj  '"e'^^j^  ' njj 


(38) 


Define  the  inner  product  of  two  vector  functions  on  S to  be  the 
Integral  over  S of  the  dot  product  of  these  two  vector  functions.  Because 
the  field  operators  in  (37)  and  (38)  are  the  same  as  those  considered  in  19], 
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only  the  nth  term  of  the  sum  (37)  or  (38)  contributes  to  the  Inne  product 
of  (37)  or  (38)  with  either  or  Hence,  the  Inner  produ  t of  (37) 

with  J*"  , 1=1, 2,... N,  and  , 1“1,2,...N,  successively  and  the  Inner  pro- 

duct  of  (38)  with  J_^^,  1“1,2,...N,  and  J_^^,  1=1,2, .,.N,  successively  gives 
the  matrix  equation 


nd-* 

(Y^  + 

ne 

aY^J) 

nd 

.tt 

nd-’ 

(-Y*^  - 
n€ 

<d> 

- 

, nd 
d 

(2^ 

ne 

"d 

- z“) 

d 

(-Z^^ 

ne 

6ne 

' "d 

1,  +1, 

+2 , . . , 

ne  n nd  ne  n nd  n 

e 'e 


for  n = 0,  +1,  +2, In  (39),  V*",  I^,  and  are  co]  imn 

— — n n n n 

vectors  of  he  coefflclen  s appearing  in  (33)  and  (34).  Also, 


nf  Lj 


. • n X )ds 

ni  “»f  Vn  1 


(zP'’) 

^ nf^•j 


ini  • 


vP  = — ff  jP  • ds 

nl  T]  ^ J J 

S 

iP,  - jP  , • n X ds 

ni  J J /“-ni 


where  p may  be  either  t or  ((>,  q may  be  either  t or  (^,  and  f may  be 
either  e or  d.  If  p-q  in  (40),  it  matters  whether  the  magnetic  field 
Hj;(^j)  is  evaluated  just  outside  or  just  inside  S.  The  y's  without 
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carets  In  (39)  are  given  by  the  right-hand  side  of  (40)  in  which  the  mag- 
netic field  is  evaluated  just  inside  S.  The  Y's  with  carets  in  (39)  are 
given  by  the  right-hand  jide  of  (40)  with  magnetic  field  evaluation  just 
outside  S. 

The  Y and  Z submatrices  on  the  left-hand  side  of  (39)  are  the  same 

as  in  equation  (88)  on  page  24  of  [9]  with  the  reservations  that  the  caret 

on  Y denotes  magnetic  field  evaluation  just  outside  S,  and  the  extra  sub 

script  e or  d denotes  radiation  in  either  y , e oi  y,,  t The  I colum 

e e d d 

vectors  on  the  right-hand  side  of  (39)  are  the  sain  as  in  equation  (88) 
on  page  24  of  [9]  whereas  the  V column  vectors  in  (39)  are  the  same  as  tie 
Vs  without  carets  in  equation  (88)  on  page  24  of  |9]. 

The  solution  and  to  the  matrix  equation  (39)  det  rmines 

n n n n 

the  equivalent  electric  and  magnetic  currents  J and  M according  to  (33  and 
(34),  From  Fig.  2,  these  currents  radiate  in  y^,  to  produce  the  scattered 
field  outside  S. 


IV.  FAR  FIELD  MEASUREMENT  AND  PLANE  WAVE  EXCITATION 

In  this  section,  meas  irement  vectors  are  used  to  obtain  the  far  field 

of  the  equivalent  surface  cu  -rents  J and  M radiating  in  y , e . Th  s far 

e e 

field  is  the  far  field  scattered  by  the  homogeneous  body  of  revolut  on.  For 
plane  wave  excitation,  the  composite  vector  on  the  right-hand  side  jf  (39) 

Is  expressed  in  terms  of  these  measurement  vectors. 

By  reciprocity, 

’ 

£"•!£=  (J(r)  • E(Iil^)  - M(r)  • H(IJl^))ds  (44) 

s 

g 

where  E is  the  far  electric  field  due  to  J and  M,  I£.  is  a receiving  elec- 
trie  dipole  at  the  far  field  measurement  point,  E(I)l  ) is  the  electric  flel  i 
due  to  IJl_,  and  H(Iil  ) is  the  magnetic  field  due  to  li  . Both  E(Ii  ) and 
H(Ii^)  are  evaluated  at  point  ^ on  S where  r is  the  point  at  which  the  dif 
ferentlal  portion  of  surface  ds  is  located.  If  is  tangent  to  the  radia  ion 
sphere. 
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/ 


*~r 


H(lO 
^ —r 


-1k.n  e 


-jkr 


rJ.  e 


Airr 

r 

-jkr 


4Tir 


(k  X 


li.  )e 

-MT 


-itr 

•HT 


(45) 


(46. 


where  is  the  distance  between  the  measurement  point  and  the  origin  in 

the  vicinity  of  S.  Also,  k^  is  the  propagation  vector  of  the  plane  wave 

coming  from  I£^,  k is  the  propagation  constant  and  n is  the  Intrinsic 

Impedance  of  the  medium  outside  S.  To  simplify  the  notation  in  this  se  - 

tion,  we  have  omitted  the  subscript  e from  all  parameters  dependent  on  i he 

medium.  It  is  understood  that  all  far  field  measurement  vectors  and  plane 

wave  excitation  vectors  depend  only  on  the  external  medium  y , e . 

e e 

Substitution  of  (33),  (34),  (45),  and  (46)  into  (44)  gives 


-jkr 


s _ -jn  e 

Ee  _ 


r ” 


“ r»  r»  n r»  r» 


jn* 


4Trr  ‘-  'nn  nn  nn.  nn' 

r n=“<» 


(47 


for  li  = u^  and 
~,r  '*•6 


-jkr 


jn*. 


f = Z2IE I _ -*e-*  ^ R**f*)^-^  r 

* 4TTr  nn  nn  nn  nn 

r n=-“ 


(48) 


r T r 

for  li  = u.  where  u„  and  u^  are  unit  vectors  in  the  0 and  * direct  ons 
~-r  «~*  r r 

respectively.  As  shown  in  Fig.  7,  0^  and  *^  are  the  angular  coordin-.tes  of 
the  receiver  location  at  which  is  placed.  In  (47)  and  (48),  E®  ind  E® 

0 9 

are  the  6 and  * components  of  E^.  Also,  V*^,  and  are  c 'lumn 
r r «-  n n n n 

vectors  of  the  coefficients  appearing  in  (33)  and  (34).  Furthermo’e,  s 

n 

a row  vector  whose  jth  element  is  given  by 


li 


li 


rP^  = k 
nj 


-jn*^  f 


4 


(•9) 


where  p may  be  either  t or  * and  q may  be  either  0 or  *.  In  view  of  (35) 
and  (36),  (49)  is  the  same  as  equation  (92)  on  page  26  of  [9].  It  is 
shown  in  [9]  that  the  right-hand  side  of  (49)  does  not  depend  on  *^. 
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For  plane  wave  incidence  and  expansion  functions  J . and  i iven 

— nj  -nj 

by  (35)  aid  (36),  the  equivalent  currents  (33)  and  (34)  and  the  fields  (47) 
and  (48)  have  special  forms.  To  obtain  these  forms,  assume  that  the  Incident 
electromagnetic  field  E^,  is  either  a 0 polarized  field  defined  by 

“jk  • r 

= kn  lu  e ^ (50) 

i -ih-  ' ^ 

H = -kUye  (51) 

or  a 4i  polarized  field  defined  by 

. -jk  • r 

E^  = kn  u e ""  (52) 

- ^y 


where  k^  is  the  propagation  vector  and,  as  shown  in  Fig.  7,  u^  and  are 
unit  vectors  in  the  6^  and  y directions  respectively.  Here,  6^  is  i he 
colatitude  of  the  direction  from  which  the  incident  wave  comes,  k^  is  in 
the  xz  plane.  No  generality  is  lost  b/  putting  ^ in  the  xz  plane  because 
if  ^ wert  shifted  out  of  the  xz  plane  by  an  angle  ^ , the  response  would 
also  be  shifted  by  the  same  angle  (ji^. 

Substituting  (50)  and  (51)  into  (42)  and  (43),  then  substituting  (52) 
and  (53)  into  (42)  and  (43),  next  taking  advantage  of  the  relationships 


J , X n = . 

ni  — nl 


j'*’  X ,1  = - 


which  are  apparent  from  (35),  (36)  and  Fig.  6,  then  comparing  the  results 
with  (49),  and  finally  using  equation  (104)  on  page  29  of  [9],  we  obtain 


1 


*t9 

V 

n 

n 

n 

-P* 

n 

n 

n 

n 

n 

|te 

n 

n 

-P^ 

n 

-^e 

n 

n 

n 

-p^ 

n 

-P^ 

n 

(56) 


:*■  if 

The  first  superscript  on  and  in  (56)  is  the  superscript  which  aDpears 


it 

on  the  right-hand  side  of  (39).  The  second  superscript  on  V and  I Ln  (56) 

n n 

denotes  the  polarization  of  the  incident  plane  wave.  If  this  second  s.uper- 
.script  is  ),  the  9 polarized  field  given  by  (50)  and  (51)  is  incident  If 
this  seconi’  subscript  is  ij),  the  (}>  polarized  field  given  ty  (52)  and  ( >3)  is 
incident.  The  jth  element  of  the  column  vector  on  the  right-hand  side 


of  (56)  is  given  by  (49)  with  6^  replaced  by  9^..  Conceding  that  9^  does  not 


appear  explicitly  in  (49) , we  really  mean  that  9^  is  replaced  by  9^.  after  the 
surface  integral  in  (49)  is  evaluated.  In  other  words,  9^  is  replaced  by  9^ 
in  equation  (95)  on  page  27  of  [9]. 


Fir  plane  wave  incidence,  the  +n  and  -n  terms  in  formulas  (33)  and 
(34)  for  the  equivalent  currents  can  be  combined  as  follows.  According  to 
equations  (102)  and  (103)  of  [9],  the  Y and  Z submatrices  in  (39)  are  either 
even  or  odd  in  n.  The  even-odd  properties  in  n of  the  submatrice;  of  the 
square  matrix  on  the  left-hand  side  of  (39)  are  tabulated  as 


where  + denotes  an  even  submatrix  and  - denotes  an  odd  submatrix.  It 
I ollows  that  the  submat  ‘ices  of  the  inverse  of  the  square  matrix  om  the 
j eft-hand  side  of  (39)  have  even-odd  properties  in  n given  by 
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J 


From  (56)  and  the  even-odd  properties  of  given  by  equation  (104)  on 
page  29  of  [9],  the  column  vectors  and  on  the  right-hand  side  of 
(39)  are  either  even  or  odd  in  n.  The  even-odd  properties  of  the  sub- 
matrices on  the  left-hand  side  of  (56)  are  tabulated  as 


Because  of  the  above  even-odd  properties  of  the  square  matrix  on  the  left- 
hand  side  of  (39)  and  the  column  vector  on  the  right-hand  side  of  (39), 
the  solutions  to  (39)  satisfy 


|te 

n n 


The  first  superscript  on  the  column  vectors  and  in  (57)  is  that 

_ ZT' 

which  appears  on  the  column  vectors  ^ and  t on  the  left-hand  side  of 

n n 

(39).  The  second  superscript  on  the  column  vectors  in  (>7)  denotes  ei  her 
the  6 or  the  ()>  polarized  incident  plane  wave.  Substitution  of  (57),  ( 5), 
and  (36)  into  (33)  and  (34)  yields 
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J*  = + I 

° n=l 

i M***  = (f^'*’)u  + I 
n o 

n=l 

for  the  (p  polarized  incident  vave.  In  (58)-(61),  f is  a row  vector  of  | 

the  f.(t).  The  superscript  6 or  <J)  on  J or  M in  (58)-(6])  differentiates  | 

^ I 

the  equivalent  currents  for  the  0 polarized  incident  wave  from  those  fc  i 

the  (|i  polarized  incident  wave.  ' 


2j(ft'^‘^)u  sin(n<t,)  + 2(ft'*’**’)u^  cos(n()))  (60) 

n n 

2(f^*^)  i cos(n<|!)  + 2j  (fV*^*^)!!  sin(n()))  (61 


The  far  scattered  field  i (47)  and  (48)  are  specialized  to  the  f 
polarized  incident  plane  wave  y appending  the  additional  subscript  6 to 

g 

E on  the  left-hand  sides  of  ( 7)  and  (48)  and  the  additional  superscript 
0 to  V^,  and  on  the  right-hand  sides  of  (47)  and  (48) . More- 

over, in  view  of  equation  (104)  on  page  29  of  [9]  and  (57),  the  +n  and  -n 
terms  in  (47)  and  (48)  can  be  combined.  As  a result,  (47)  and  (48)  become 


4TTr 

r 


00 


n n n n 


+ R'<>®t'»>®) 

n n n n 


cos(n(|»^)  j- 


2'rTr 

r 


I R 

n n ’ 

n=l 


^t4^te  ^ 
n n 


(62) 


+ R^‘*^*^)8ln(n(t>^)  (63) 

for  the  0 polarized  Incident  plane  wave.  Similarly,  (47)  and  (48) 
become 
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(64) 


-jkr^ 


=.s  _ ne 


2irr 


r n=l 
jkr 


I + R^^**)sin(n<t.  ) 

nn  nn  nn  nn  r 


E®  = Ijns 1 \ - R^^^^  + + 2 I (-R^^ 

dxli  Airr  Loo  oo  nn 

r n=l 


n n 


~t^t*  + R*'^t^‘'>)cos(n4.  )]  (63) 

n n n n r J 


g 

for  the  (|>  polarized  incident  plan'  wave.  The  first  subscript  on  E on 
the  left-hand  sides  of  (62)-(65)  denotes  the  receiver  polarization  and 

g 

the  second  subscript  on  E denotes  the  transmitter  polariz  ition. 

The  scattering  cross  section  a is  defined  by 

pq 

4Trr^|E® 


pq 


,i|2 


(66) 


where  p is  e ■ ther  0 or  cji  and  q is  either  6 or  cji.  In  (66),  E^^  is  a com- 
ponent of  tht  scattered  field  given  by  (62)-( j5)  and  is  the  magnitude 

of  the  electric  field  of  the  inci  lent  plane  wave.  According  to  (50)  and  (52), 


for  both  polarizations  so  that 


= kn 


Airr^lE® 


r'  pq' 


pq 


1,2  2 

k n 


Normalized  versions  of  (68)  are 


a 4r^ I E®  I ^ 

pq  = r'  pq' 

2 ,222 

ira  k a n 


(67) 


(68) 


(69) 


2|„s  I 2 
a r E 
pq  _ r'  pq' 

.2  ~ 2 
A irri 


(70) 


where  a is  some  characteristic  length  associated  with  the  scat ter or  and  X is 
the  wavelength  in  the  external  medium. 
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V.  EXAMPLES 


A computer  program  has  been  written  to  calculate  the  equivalent 
currents  and  scattering  patterns  for  a dielectric  body  of  revolution  ex- 
cited by  an  axially  incident  plane  wave.  This  progran  Is  described  and 
listed  in  Part  Two.  Some  computational  results  obtained  with  this  program 
are  given  in  this  section. 

Figures  8 and  9 show  the  ni;>gnltude  and  phase  of  the  normalized 
Jq  J.  Mq  m 

equivalent  currents  — , and  ^ on  the  surface  of  a dielectric 

n H L 

y y X X 

sphere  for  which  ka  = 3 and  = 4.  Here,  k is  the  propagation  constan 

in  free  space,  a is  the  radius  of  the  sphere  and  is  the  relative 

dielectric  constant  of  the  sphere.  Figure  8 represents  our  solution  o the 
PMCHW  formulation.  Figure  9 depicts  our  solution  of  the  Muller  formulation. 


In  Figs.  8 and  9,  the  incident  field  is  a plane  wave  traveling  in 

the  positive  z direction.  THETA  = 0“  is  the  forward  scattering  direction 

and  THETA  = 180°  is  the  backscatte i ing  direction.  The  incident  field  is 

given  by  (50)  and  (51)  with  6 = 180°.  The  origin  r = 0 is  at  the  center 

t ^ 

of  the  sjthere.  In  F gs.  8 and  9,  J Is  the  u = - u component  of  electric 

current  (58)  versus  ' in  the  (p  = 0 plane,  J,  Is  the  u.  component  of  (58) 

versus  0 in  the  p = 90°  plane,  M.  is  the  u„  = - u^  component  of  ma^^etlc 

current  (59)  versus  0 in  the  i(i  = 90°  plane,  and  M is  the  u,  component  of 

(59)  versus  0 in  the  4)  = 0°  plane.  For  axial  incidence,  only  the  n=l  term 

is  present  in  (58)  and  (59) . The  symbols  x and  + denote  respectively 

magnitude  and  phase  of  the  method  of  moments  solution  for  the  pertinent 

component  of  the  electric  or  magnetic  current.  The  solid  curves  . re  the 

exact  equivalent  currents  obtained  from  the  Mle  series  tolutlon  [11].  The 

normalizing  constants  E and  H are  defined  in  terms  of  the  Incident  field 
*=*  X y 

(50)  and  (51)  by 


E 

X 


u 

''X 


H 


y 


u 


'-y 


(71) 


[11]  R.  F.  Harrington,  Time-Harmonic  Electromagnetic  Fields,  McGrtw-Hill 
Book  Co.,  1961.  Section  6-9. 
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where  k and  n are  respectively  the  jropagatlon  constant  of  fre«  space  and 
the  Impedance  of  free  space. 

The  currents  of  Figs.  8 and  9 were  obtained  by  using  a 20  point 
Gaussian  quadrature  formula  for  all  integrations  in  ip.  All  integrations 
over  the  functions  {fj(t:)}  in  t • aCtt-B)  were  done  by  sampling  each 
fj(t)  four  times.  The  fj(t)}  consisted  of  14  overlapping  triangle 
(divided  by  the  cylindrical  coordinate  radius)  functions  equally  spaced 
in  6.  More  precisely, 

NP  = 31 

NPHI  = 20  (72) 

MT  = 2 

where  the  above  variables  are  input  data  for  the  computer  program  described 
and  listed  in  Part  Two,  Section  V. 

Figures  10  and  11  show  the  scattering  patterns  radiated  by  the  cur- 
rents of  Figs.  8 and  9 respectively.  The  symbols  x and  + denote 

^60  ^ 1>0 

— 2 and  — ^ respectively.  The  solid  curves  are  the  exact  patterns  obtained 
wa  Tra 

from  the  Mle  series  solution  [11].  The  patterns  o.^  and  a,,  are  given  by 

00  0 

(69),  (62)  and  (63).  Here  is  the  0 polarized  pattern  versus  0^  in  the 
41  = 0 plane  and  is  the  4 polarized  pattern  versus  0^  in  the  = 90° 
plane.  The  THETA  in  Figs.  10  and  11  refers  to  0^.  For  axial  incidence, 
only  the  n=l  terms  are  present  in  (62)  and  (63).  Elsewhere  [3,  12],  the 
pattern  Ogg  is  called  the  horizontal  polarization  because  it  is  polarized 
parallel  to  the  scattering  plane.  Similarly,  the  pattern  a.,  is  called  the 
vertical  polarization  because  it  is  polarized  perpendicular  to  the  scatter- 
ing plane. 

Figures  12-17  show  the  scattering  patterns  for  three  other  d: electric 
spheres.  Figures  12  and  13  are  for  relative  dielectric  constant  = 1.1, 
Figs.  14  and  15  for  e « 10.,  and  Flgj.  16  and  17  for  e = 20.  All  other 

r o j. 

[12]  P,  Barber  and  C.  Yeh,  "Scattering  of  Electromagnetic  Waves  by 

Arbitrarily  Shaped  Dielectric  Bodies,"  Applied  Optics,  vol.  14,  No  12, 
December  1975,  pp.  2864-2672. 


iiation  respectively.  Solid  line  denotes  exact  solu  Ion. 
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Mlilier  solution.  Symbols  ' auJ  + denote  horizontal  polai  ization  and 
vertical  polarization  respectively.  Solid  line  denotes  exact  solution. 


27 


O.OOOH 1 1 m I 1 

0 30  60  90  120  150  1 


THETfl 

Kip,.  12.  Plane  wave  scattering  p.  tterns  for  dielectric  sphere,  ka  “ 3,  ‘ = 1 . 

PMCHW  solution.  Syndsolt^  * and  + denote  horizontal  pol  irizati  m'^and 
Vertical  polarization  respectively.  Solid  line  denoteb  exact  solution 
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fig-  13.  Plane  wave  scattering  patterns  for  dielectric  sphere,  ka  = 3,  =11 

Muller  solution.  Symbols  x and  + denote  horizontal  polarization  and 
vertical  polarization  tespei  tlvely . Solid  line  denotes  exact  aol ui i on 


(•’ig.  16,  Plant'  wave  scattering  patterns  for  dielectric  sphere,  ka  = 3,  c = 2(  , 
PMCIIW  solution.  Symbt  Is  x and  + denote  horizontal  polarization  and 
vertical  polarization  i espectlvely . Solid  line  denotes  exac  solui i 'n 


SIGMfl/iTR2 


100 


parameters  in  Figs.  12-17  are  the  same  as  in  Figs.  10  and  11.  In  Figs.  12 
and  13,  values  less  than  0.0001  are  plotted  at  0.0001. 

Figures  18  and  19  show  the  computed  scattering  patterns  of  a finite 
dielectric  cylinder  of  radius  a and  height  2a  when  a is  0.25  free  space 
wavelengths.  The  relative  dielectric  constant  of  the  cylinder  is  i ^ 

The  incident  field  is  a plane  wave  traveling  in  the  positive  ?.  direction, 
the  same  field  which  was  incident  upon  the  previous  dielectric  spheres. 

^66  6 

Figure  18  shows  — ^ and  — ^ as  obtained  from  our  solution  ol  the  PMCHW 
ira  Tra 

Oge  a 

formulation.  Figure  19  shows  — j “2  obtained  from  our  solution 

tra  Tra 

^09  ^(bS 

of  the  Muller  formulation.  The  patterns  — ^ and  are  plotted  with  the 

ira  ira 

symbols  x and  + respectively. 

The  equivalent  currents  which  radiate  the  patterns  of  Figs.  18  and  19 
were  obtained  by  using  a 48  point  Gaussian  quadrature  formula  for  all  inte- 
grations in  (j>.  All  integrations  in  t over  the  functions  {fj(t)}  were  done  by 
sampling  each  four  times.  The  {f^(t)}  consisted  of  11  overlapping 

triangle  (divided  by  the  cylindrical  coordinate  radius)  functions  equally 


as  obtained  from  our  solution 


spaced  in  t.  More  precisely. 


Nl’  = 25 


NPHl  = 48 


MT  = 2 


where  the  above  variables  are  input  data  for  the  computer  program  described 
and  listed  in  Part  Two,  Section  V. 


VI.  DISCUSSION 


According  to  Figs.  12  and  13,  the  scattering  patterns  obtained  from 
our  solution  of  the  Mliller  formulation  are  more  accurate  thai  those  ob- 
tained from  our  solution  of  the  PMCHW  formulation  for  the  di« lectrlc  sphere 
with  ka  = 3 and  = 1.1.  From  plots  not  Included  in  this  n port,  we  ob- 
served that  both  our  PMCHW  solution  and  our  Mliller  solution  Cor  the  equivalent 
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currents  on  the  dielectric  spl  ere  were  reasonably  accurate.  However,  tie 

following  argument  shows  that  when  is  near  one,  a slight  inaccurac}  in 

the  equivalent  currents  could  affect  the  scattering  patterns  drastically. 

As  e approaches  one,  the  equivalent  electric  and  magnetic  currents 
^ i i 

approach  n x H and  E x n respectively  whereas  the  scatteming  patterns 
approach  zero.  This  means  that  the  equivalent  electric  and  magnetic 
currents  produce  fields  which  nearly  cancel  each  other.  Hence,  a slight 
inaccuracy  in  the  equivalent  currents  coul(  cause  a large  percentage 
inaccuracy  in  the  scattering  patterns. 

We  believe  that  our  Muller  solution  is  more  acc  irate  than  our 
PMCHW  solution  whenever  is  close  io  one.  When  a and  B are  given  by 
(22)  and  (23)  as  in  the  Muller  formulation,  the  left-hand  sides  of  (9) 
and  (10)  approach  -M  and  ^ respective  y as  approaches  one.  In  this 
case,  the  expected  solution 

J = n X 
M = E^  X n 

can  b.‘,  obtained  by  Inspection  of  (9)  and  (10).  However,  if  a = 6 = L 
as  in  the  PMCHW  foirmulation,  the  solution  to  (9)  and  (10  is  not  obvious 
when  = 1 because  the  field  operators  on  the  left-han  I sides  of  (9) 

and  (10)  are  not  diagonal.  With  our  Muller  solution,  tl  ^ matrix  on  the 
left-hand  side  of  (39)  would  become  tridiagonal  for  1 if  its  first 

two  rows  of  submatrices  were  interchanged.  With  our  PM  HW  solution,  no 
such  simplification  of  this  matrix  is  posslbl j for  = 1. 

We  recommend  at  least  10  expansion  functions  per  wavelength  per 
component  along  the  generating  curve  of  the  dielectric  body  of  revolution. 
For  example,  if  the  generating  curve  were  one  wavelength  long,  the  order 
of  the  square  matrix  on  the  left-hand  side  of  (39)  should  be  at  least  36. 
The  number  36  is  arrived  at  as  follows.  There  should  be  at  least  9 expan- 
sion fund  Ions  per  component  of  current.  We  say  9 expansion  functions 
rather  than  10  because  we  are  using  overlapping  triangle  function!  with  no 
peak  of  triangle  function  at  either  ends  of  the  generating  curve.  There 
are  two  components  of  electric  current  and  two  components  of  magnetic 
current. 
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According  to  equations  (20) -(2 3)  of  [9]  and  (58) -(61)  of  [91,  each 
element  of  the  square  matrix  in  (39)  is  a triple  integral  consist!  ig  of 
one  integration  with  respect  to  ij>  and  two  integrations  with  respect  to  t. 
The  (()  integral  is  evaluated  by  using  a Gaussian  quadrature  formula.  Each 
t integration  is  done  by  crude  sampling  akin  to  the  trapezoid  rule.  In 
any  case,  there  should  be  at  least  10  samjle  points  per  wavelength  in  the 


media  in  question.  For  instance,  if  is  the  largest  cylindrical 

coordinate  radius  of  the  dielectric  body  of  revolution  and  t is  the 

max 


length  of  the  generating  curve,  t'len 


where  NP  and  NPHI  are  input  arguments  of  the  subroutine  YZ  described  md 
listed  in  Part  Two,  Section  II  and  X is  the  wavelength  in  the  media  n 
question.  If  f = e in  (40)-(41)  then  X is  the  wa/elength  in  the  exi  jmal 
media,  but  if  f = d in  (40) -(  tl)  then  X is  the  wavelength  inside  th  ; dif- 
fracting body  of  revolution.  The  main  program  in  Part  Two,  Section  V is 
oversimplified  in  that  it  uses  the  same  values  of  NP  and  NPHI  for  both 
f = e and  f = d. 


Loss  of  accuracy  in  the  computed  patterns  of  Figs.  16  and  17  may 
be  due  to  the  fact  that  (74)  was  violated.  According  to  (74),  the  values 
of  NP  and  NPHI  for  f = d should  be  nearly  70  or  greater  Insteac  of  the  low 
values  appearing  in  (72).  Unfortunately,  increasing  the  valuet  of  the 
variables  NP  and  NPHI  Increases  the  computer  time  required  to  solve  the 
problem. 


We  have  been  trying  to  obtain  accurate  numerical  results  for  the 
dielectric  sphere  for  which  a = 0.2  free  space  wavelengths  and  = 80. 
from  our  g.  neral  dielectric  body  of  revolution  program.  We  have  n >t  been 
able  to  obtain  clear-cut  convergence  with  respect  to  the  variables  on  the 
left-hand  sides  of  (72)  because  we  could  not  afford  to  Increase  them  as 
much  as  desired.  Our  PMCHW  solution  and  our  Muller  solution  for  the 
equivalent  currents  and  scattering  patterns  differ  from  each  other  and 
from  the  exact  solution. 
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For  the  sphere,  each  element  of  the  quare  matrix  on  the  left-hand 
side  of  (39)  can  be  written  as  a sum  >ver  the  Infinite  set  of  spherical 
inodes.  So  far,  we  have  not  been  able  to  successfully  Implement  thi i altet  - 
nate  evaluation  of  the  matrix  elements  In  terms  of  spherical  modes.  The 
major  difficulty  seems  to  be  lack  of  agreement  of  a few  matr  x elements  for 
which  both  expansion  and  testing  fun  :tlons  are  near  one  of  the  poles  of  the 
sphere. 

Both  the  PMCHW  solution  and  the  Muller  solution  are  obtained  by  taking 
a linear  combination  of  (5)  and  (7)  and  a linear  comblnatlo i of  (6)  and  (8  . 
There  are  two  other  possibilities  which  are 

(1)  A linear  combinat  Lon  of  (5)  and  (6)  and  a linear  combi  - 
nation  of  (7)  am:  (8). 

(2)  A linear  combination  of  (5)  and  (8)  and  a linear  combi- 
nation of  (6)  and  (7) . 

These  other  two  possibilities  give  rise  to  alternative  numerical  solutioi  s 
which  may  compare  favorably  with  the  PMCHW  solution  and  the  Miiller  solu  ion. 


APPENDIX  A 


THE  EQUIVALENCE  PRINCIPLE 


Let  be  an  electromagn^stic  f eld  defined  outside  a closed 

surface  S.  The  perroi  ability,  pennlttlvlt> , and  electric  and  magnetic 
source  currents  outside  S are  and  respectively.  Thi  ■ 

outside  situation  where  the  subscrlf  test  ands  for  "exterior  medlui  " 

Is  shown  in  Fig.  A-1.  In  Fig.  A-1,  the  media  and  sources  Inside  S are 
undlsclosei' . Let  E^,  be  an  electromagnetic  field  defined  Inslda  S 
where  the  p ;rmeablllty,  permittivity,  and  electric  and  m ignetlc  sources 
are  and  respectively.  This  'nslde  situation  where  the 

subscript  d stands  for  "diffracting  mi>dium"  is  shown  in  Fig.  A-2.  In 
Fig.  A-2,  the  media  and  sources  outsice  S are  undisclosed.  The  equi- 
valence principle  states  that  the  solution  to  the  composite  radiation 


problem  consisting  of  medium  y , e an<l  sources  J , M outside  S,  medium 

e e ^ ^ 

y^  , and  sources  Inside  S,  am  electric  and  magnetic  surface 

currents  J[,  M on  S given  by 

J = n X (H  - 

(A-1)  1 

M = (E  - E_.)  X n 
^ ^•*<1 

j 

(A-2)  i 

1 

where  n is  the  exterior  unit  normal  vector  on  S is  the  ccmposlte  elec- 
tromagnetic field  E,  H defined  by 

1 

1 

E,  H = E , H outside  S 

(A-3) 

E,  H = Ej,  H , Inside  :! 

(A-4) 

The  composite  radiation  problem  is  sh  )wn  In  Fig.  A-3  which  Is  entitled 

composite  situation. 

Tiie  equivalence  principle  is  proved  by  showing  that  the  configu- 
ration of  media  and  sources  in  Fig.  A-3  gl\es  rise  to  the  composite  field 

E,  H defined  by  (A-3)  and  (A-4) . Now,  E,  H will  be  the  field  gem  rated 
by  the  media  and  sources  of  Fig.  A-3  if  E,  H satisfies  Maxwell's  t quatlons 
with  source  terms  Included  and  the  radlat  on  condition  at  Infinity.  E,  H 


satisfy  Maxwell's  equations  ou  side  S and  inside  J respectively. 
E,  H also  satisfies  the  radiation  conlition  at  infinity  because  E , H 
satisfies  the  radiation  condition  at  infinity.  It  remains  to  be  shown 
that  Maxwell's  equations  for  E,  H exhibit  the  surface  current  sources 
J and  M given  by  (A-1)  and  (A-2) . 

It  is  well  known  that  a surf  ice  current  source  on  S gives  rise  o 
a discontinuity  across  S of  t le  tangential  component  of  the  field.  The 
prect  ding  statement  is  easily  verified  by  means  of  an  argument  based  m 
the  iategral  forms  of  ^axwell's  equations.  Now,  this  sam ^ argument  < in 
be  CO  istrued  to  imply  that  a discontinuity  across  S of  the  tangentla 
compo  ent  of  the  field  gives  rise  to  a surface  current  source  on  S. 

Hence  Maxwell's  equations  for  E,  exhibit  the  electric  and  magneti' 
surface  current  sources  ^ and  M^on  S given  by  (A-1)  and  (A-2) . There- 
fore, E,  H is  the  solution  to  the  couposlte  radiation  problem  shown  in 
Fig.  A-3  because  E,  H satisfies  Maxw-ll's  equations  with  sources  and  the 
radiation  condition  at  infinity. 


APPENDIX  B 


PROOF  THAT  THE  SOLUTION  TO  E(  UATIONS  (5) -(8)  IS  UNIQUE 


The  solution  J,  M to  (5)- (8)  will  be  unique  If  the  assoclatt d set 
of  homogeneous  equations 


-n  X E =0  (B-1) 
x-e 

-n  X H~  = 0 (B-2) 
-n  X st  = 0 (B-3) 
-n  X nt  = 0 (B-4) 


has  only  the  trivial  solution  J = M = 0. 

From  (B-1)  and  (B-2),  J,  M radiate  In  to  produce  a field 

whose  tangential  components  are  zero  just  Inside  S.  Hence,  according 

to  the  re 1 It Ion  between  J,  M and  the  discontinuity  of  tangential  field 

across  S a;  exemplified  by  (A-1)  and  (A-2),  the  field  ^ radiated  by 

J,  M in  p , E outside  S satisfies 
e e 


n X H 


(B-5) 


E X n = M 
•xe  — 


(B-6) 


just  outside  S.  See  Fig.  B-1. 


From  (B-3)  and  (B-4)  the  electric  and  magnetic  currents  -J,  -M  radiate 

In  p,,  E , to  produce  a field  whose  tangential  components  are  zero  just  outside 
d d 

S.  Hence,  according  to  the  relation  between  -J,  -M  and  the  discontinuity  of 

tangential  field  across  S,  the  field  -E,,  -H.  radiated  by  -J,  -M  In  p.,  e. 

#^d  ^ d d 

satisfies 


n X (-Hj)  = J 

(B-7) 

(-E.)  X n = M 

4^d  ^ ^ 

(B-8) 

just  Inside  S.  See  Fig.  B-2. 

I'he  equivalence  principle  Is  used  to  combine  the  outside  situation 
In  Fig.  }-l  with  the  inside  situation  In  Fig.  B-2  to  obtain  the  composite 
sltuatioi  shown  in  Fig.  B-3.  Because  of  (B-5)-(B-8),  the  compt site 
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situation  In  Fig.  B-3  is  source -free.  Therefore,  the  field  in  Fig.  B-3 
is  zero  everywhere.  Hence,  the  fields  in  Figs.  B-1  and  B-2  ire  zero 
everywhere  which  implies  that  J = M = 0. 

Thus,  the  solution  to  (5)-('8)  is  unique  because  the  associate  1 set 
(B-1)-(B-4)  of  homogeneous  equations  has  only  the  trivial  solution. 
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APPEND  X C 
MORE  EX\MPLES 

The  equivalent  currents  and  scattering  patterns  for  the  dielectric 
spheres  for  which  = 4 and  ka  = 4,  5,  and  6 are  plotted  in  Appendix  C. 
Figures  C-1  to  C-4  are  for  ka  = 4,  Figs.  C-5  to  C-8  are  for  ka  = 5,  and 
Figs.  C-9  to  C-12  are  for  ka  = >.  All  other  parameters  are  the  same  as  in 
Figs.  8-11  in  Section  V.  In  particular,  the  input  data  for  the  comp  iter 
program  which  generated  the  method  of  moments  results  plotted  in  Figs.  C-1 
to  C-12  is  given  by  (72) . 

It  is  evident  from  Figs.  8-11  and  Figs.  C-1  to  C-12  that  the  method 
of  moments  solutions  for  the  equivalent  currents  and  scattering  patte  ns 
are  not  as  accurate  at  ka  = 4,  5,  and  6 as  at  ka  = 3.  Los i of  accura  y 
at  the  higher  values  of  ka  may  be  due  to  violation  of  (73) . 
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Fig.  C-3.  Plane  wave  scattering  patterns  for  dielectric  sphere,  ka  = A,  e 4,  PMCHW 
solution.  Symbols  x and  + denote  horizontal  polarization  and  vert  Leal  polari 
zation  respectively.  Solid  line  denotes  exact  solution. 
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l ig.  C-4.  I’lai  ■>  wave  scattering  patterns  for  rlleleci  ric  sphere,  ka  * e 4, 
MiilJer  solution.  Symbols  and  + denote  horizontal  polarlza  lon'^arJ 
vertical  polarization  respectively.  Solid  line  denotes  exac  solution. 


50 


180 


IGMH/nflS 


f--7.  Plane  wave  scattering  patterns  For  dielectric  sphere,  ka  = 5,  r - 4,  PMCHW 
solution.  Symbols  ''  and  + denote  horizontal  polarization  md  ve’^tical 
polarization  respectively.  Solid  Hie  denotes  exact  solu  ion. 
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C-9.  Equivalent  electric  and  magnetic  currents  for  dielectric  sphere,  ka  = 6,  tj.  “ PMCHW  solution 
Symbols  x and  + denote  magnitude  and  phase  respectively.  Solid  line  denotes  exact  solution. 


Fig.  C 11.  F lane  wave  sc.ittering  patterns  for  dirlectiic  sphere,  ka  - 6, 

PMCHW  solution.  Symbols  ^ ami  + denot«  horizontal  polarization  and  -'erti 
cal  polarization  respectively.  Solid  line  denotes  exact  solution. 
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1 . INTRODUCTION 

The  computer  program  calculates  the  equivalent  electric  and  magnetic 
currents  (58)  and  (59)  and  the  scattering  patterns  (70)  for  a loss-free 
homogeneous  dielectric  and/or  magnetic  material  body  of  revolution  immersed 
in  an  axially  incident  plane  wave.  Tliis  computer  progra  i consists  of  a 
main  prO};ram  and  the  subroutines  YZ,  PLANE,  DECOMP  and  S )LVE. 

Part  Two  consists  of  definitions  of  the  input  and  output  for  the 
subroutines  YZ,  PLANE,  DECOMP,  and  SOLVE,  listings  of  these  subroutines, 
definitions  of  the  input  and  output  for  the  main  program,  a verbal  flow 
char';  of  the  main  program,  and  a listing  of  the  main  program  with  sample 
input  and  output.  The  subroutines  YZ  and  PLANE  are  similar  to  subrouti  les 
of  the  same  name  in  [10].  The  subroutines  DECOMP  and  SOLVE  are,  ixcept 
for  dimension  statements,  exactly  the  same  as  in  [13].  Hence,  the  Insides 
of  the  subroutines  YZ,  PLANE,  DECOMP,  and  SOLVE  are  not  described  in  de  ail 
in  Part  Two.  Because  these  subroutines  are  quite  complicated,  a black  >ox 
approach  is  suggested  wherein  the  user  is  concerned  with  just  the  input  and 
output  of  these  subroutines.  However,  the  user  is  encouraged  to  delve 
inside  the  main  program  and  to  make  any  changes  therein  that  he  deems  neces- 
sary to  suit  his  needs. 


II.  THE  SUBROUTINE  YZ 


Description: 

The  subroutine  YZ(NN,  NP,  iPHI,  M,  MT,  RH,  ZH,  X.  A,  Y,  Z)  stores  the 

matrices  Y ^ and  Z ^ defined  by 
nf  nf 


nf 


yV 

n 


J- 


nf 


nf 


(75) 


"nf 


ni 


nf 


"nf 


Z 


(76) 


1 

I 
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by  columns  in  Y and  Z respectively.  The  submatrices  on  the  right-hand 
sides  of  (75)  and  (76)  are  given  b^  (40)  and  (41).  The  fir  it  9 argu- 
ments of  YZ  are  input  variables.  1 xcept  for  the  new  input  variables  M 
and  MT,  the  subroutine  YZ  is  the  same  as  the  old  subroutine  YZ  on  pages 
17-21  of  [10].  If  M = - 1 and  MT  = 2,  these  subroutines  are  exactly  the 
same  as  far  as  the  calculation  of  Y and  Z in  terms  of  the  rest  of  tlie 
input  variables  is  concerned. 


\ M = - 1 for  field  evaluation  just  inside  S and  M = + 1 for  field 

evaluation  just  outside  S.  M = - 1 if  f = e in  (75)  and  M = + 1 if 
f = d in  (75).  The  value  of  M is  not  used  in  calculating  (76)  because  the 
tangential  components  of  the  elJctric  field  operator  in  (41)  .\re  con- 
tinuous across  S.  All  numerical  integrations  over  t of  fj(t)  appearing 

in  (35)  and  (36)  are  done  by  sampling  each  f (t)  2*MT  times.  The  rep  e- 

d ^ 

sentations  of  pf^(t)  and  — (pf^(t))  given  by  (66)  and  (67)  of  [9]  ar ^ 


replaced  by  representations  which  contain  2*MT  Impulse  functions  ins  ead  of 
4 impuls.'  functions.  For  instance,  (66)  of  [9]  is  replaced  by 


2*MT 

pfj(t)  - iT  I T 

P=1 


p+(i-l)*2*MT 


'^^*^“%+(i-l)*MT^’ 


(77) 


i=l,2,..,N  where  N will  be  defin  :d  in  the  paragraph  which  follows  the  next 
paragrapfi.  The  T's  appearing  in  (77)  will  be  defined  by  (78). 

Seven  of  the  input  variables  are  the  same  as  in  • ne  old  subroutf  le 
YZ  on  pages  17-21  of  [10].  These  varial  les  are  defined  in  terms  of  va  iables 
appearing  in  [9]  by 

NN  = n 

NP  = P,  page  9 

NPHI  = N^,  page  13 

RH(1)  = kp^,  page  9 

ZH(i)  = .z^,  page  9 

X(k)  = page  13 

A(k)  = Aj^,  page  13 
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In  sunmary,  n denotes  the  dependence  In  (35)  and  (36),  (pT,  z ), 

are  coordinates  on  the  generating  curve,  the  k w’  ich  multi- 
plies and  is  the  propagation  constant,  and  ;.nd  Aj^  ire  th£ 
abscissas  and  weights  for  the  point  'laussian  quadrature  integration 
in  (j).  Note  that  NP-1  should  be  an  inte  ;er  multiple  of  MT.  If  NP-1  is 
not  an  integer  multiple  of  MT,  the  program  will  ignore  RH(m)  and  ZH(m) 
for 

[(NP  -1)/MT]*MT+1  < m _<  NP 

where  ((NP-1)/MT]  is  the  largest  integer  which  does  not  exceed  (NP-1)/MT. 

Minimum  allocations  are  given  by 
COMPLEX  Y(4*N*N),  Z(4*N*N) 

DIMENSION  RH(NP),  ZH(NP),  X(NPHI),  A(NPHI),  D(NG) 

PD(NG),  TP(2*MT*N),  CR(NPHI),  Cl(NPHI),  C2(NPHI), 
C3(NPHI),  C4(NPHI) 

COMMON  RS(NG),  ZS(NG),  SV(NG) , CV(NG),  T(2*MT*N) 

where 

N = [(NP-1)/MT]  - 1 
NG  = (N+1'*MT 


The  variables  in  common  make  the  results  of  some  intermediate  calculations 
done  in  YZ  available  to  the  subroutine  PLANE  described  in  Section  III  of 
Part  Two. 


We  mention  a few  portions  of  YZ  which  differ  from  the  subroutine 
listed  on  pages  18-21  of  [10].  Equation  (29)  of  [9]  has  been  generalized  to 


'^2*MT*(J-1)+I  Aj^  *W*(J-l')+q  ■ 2 ‘Sn’*(J-l)+l^‘S{T*(J-l)+I 

^ 1 ^ 

'^2*MT*J-MT+I  ° aJ  ^^2  ‘^'fr*J+q  2 ‘SiT*J+I^ ‘SlT*J+I 

for  J ” 1,2,...N  and  I = 1,.!,...M''’  where 


(78) 

(79) 


MT 

^1  ° ‘Vr*(j-i)+i 


MT 

^2  “ J-  ‘HiT*J+I 

1*1 


(80) 

(81) 


Ik 


1 
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and  the  d's  are  given  by  equation  (28)  of  [9].  Note  that  is 
electrical  length  of  generating  curve  over  which  the  first  half 
fj(t)  exists  and  that  is  the  electrical  length  of  gener;  ting 
over  which  the  second  half  of  fj(t)  exit.ts.  The  generalization 
equation  (68)  of  [9]  is 

- ^V*(J-1)-H 

2*MT*(J-1)+I  A^ 

T'  = ~**^^T*J+I 

2*MT*J-MT+I  A^ 

for  J = 1,2,..,N,  and  I = 1,2,...MT  where  A^^  and  A^  are  given  by  (80) 
and  (81) . Expressions  ( '8)  - (83)  are  calculated  in  DO  loop  68.  D' 
loop  12  accumulates  A^^  in  DEL.  DO  loop  19  puts  (78)  in  T(2*MT*(J-1) fl) 
and  (82)  in  TP(2*MT*(J-1)+I) . IX)  loop  15  accumulates  A^  in  DEL.  D<> 
loop  16  puts  (79)  in  T(2*MT*I-MT+I)  and  (83)  in  TP(2*^fr*J-MT+I) . 

Th2  subscripts  KT,  LT,  and  J1  lns'*.de  DO  loop  32  are  obtained  as 
follov.'s.  Since  the  generating  curve  consists  of  NG  = (N+l)*Mr  small 
intervals,  it  is  composed  of  (N+1)  large  intervals  where  the  mth  large 
interval  consists  of  the  ( (m-1) *MT+1) th  through  the  (m*MT)th  small 
intervals.  The  index  I of  DO  loop  60  denotes  the  Ith  small  interval. 

The  Ith  small  Interval  is  contained  in  the  (19+1) th  large  interval  where 
19  = [(I-1)/MT].  The  second  half  cf  f^g(t)  and  the  first  half  of  f^,  ^^^(t) 
are  in  this  large  interval.  The  index  K of  DO  loop  32  denotes  ]^(t)- 

Since  T((m-1)*2*MT+1)  through  T(m*2*MT)  is  allotted  to  f^(t) , in=l,2,  ..,N, 
f^(t)  is  preceded  by  (m-1)  overlaps.  For  each  overlap  the  subscript  of  T 
increases  by  an  amount  MT  not  accounted  for  by  I.  Hence,  replacirg  m by 
I9+K-1,  we  arrive  at  the  subscript 

KT  = I + (T9+K-2)*MT 

for  T.  Here,  KT  is  the  field  subscript  which  refers  to  the  testing 
function.  By  analo>y,  the  source  subscript  LT  which  refers  to  the  expan- 
sion function  is  given  by 

LT  «=  J + (,I9+L-2)*MT 


the 

of 

curve 

of 

(82) 

(83) 
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Retaining  testing  function,  we  take  the  analogously 

subscripted  function  ex  lansion  function  and  arrive 

at  the  matrix  subscript 

J1  = (J9+L-2)*N2  + I9+K-1 


where,  as  in  the  program,  N2  = 2*N. 


DO 


tt. 


loop  17  accumulates  in  R1  the  c ontribution  to  (Y  ) _ , of 

n «J ) J 


equation  (31)  of  [9]  due  to  the  ij-j— term  in  equation  (32)  of  [9  . 

This  c'ontribution  is  given  by 
2*MT 


R1 


I T 


* T, 


1 = 1 


I+2*MT*(J-1)  I+2*MT*(J-1) 


* PD(I+'1T*(J-1)) 


where,  as  in  the  program. 


PD(i)  = - M * ( 


k^d. p . 


-) 


Here,  the  factor  -M  not  included  in  enuation  (32)  of  [9]  provides  fo’" 
the  choice  of  field  evaluation  either  outside  or  inside  S. 

DO  loop  18  accumu]ates  in  R1  the  contribution  to  (¥*■*■),  , , of 

_ n J-1,J 

[31)  of  [9]  due  to  th 
This  contribution  is  given  by 


equation  (31)  of  [9]  due  to  the  term  in  equation  (31)  of  [9], 

k d^p^ 


MT 


R1  = 5;  ^I+Z*MT*(J-1)-MT  ^I+2*MT*(J-1)*  P>^(I+MT*(J-1)  ) , 

1=1 


LISTING  C:F  the  subroutine  Y7 

SUB  OUTlNF  YZ(NN,NP, NPHI ,MtMT,RH»ZHfX,AtY,Z) 

COMPLEX  U»Y (784) tZ ( 784 ) , G1 , G2, G3 f G4,G5 t G6 , Y1 , Y2, Y3, Y4 , Z 1 , Z2 , Z3 ,Z4 
DIMENSION  RH(161)  ,ZH(161),X(48),A(48)  tD(l60)  fPO(  1601  t TP(  320) 
DIMENSION  CR (48) ,C 1(48) ,C2( 48) ,C3( 48) ,C4(48) 

COMMON  kS( 160)  ,ZS( 160)  ,S V( 16  0) t CV( 160) tT(320 ) 

PI=3. 141593 
f)i  M=-M*PI 
N=(NP-1)/MT-1 
N2=2*N 
NG=( N+1 )*MT 
NGM=NG-MT 
MT2=MT*2 
D(  5 7 I = L , NG 


) 


I 
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12=1+1 

DR=FH(  12  )-RH(I) 

DZ=ZH( 12 )-ZH( I ) 

D( 1 1 =^QPT(DR*OP+DZ*OZ» 

RS(I ) .?*(RH( I2)+RH(I» ) 

ZS(1  ) .5<=(Zh(I2)+ZH(  I>  ) 

SV(I ) = nR/D( I ) 

CV(I  )=DZ/D(  I) 

PCM  I )=PIK/(D(  1 j’^RSd  ) 

57  CCNTINUE 
Jl  = 0 
J5=0 

DC.:  6 8 J = 1,N 
DEL=0. 

DC  12  1=1, MT 
J1=J1+1 
DFL=DEL+D( J1 ) 

12  continue 
ji  = n-MT 
SN=  ) . 

DC  .9  I=1,MT 
J5-  = .l5+l 
Jl=.  1+1 
SN=C  N+D( Jl ) 

TP(  J5J=D(  JD/UEL 
7(  J5)  = (Sr-.5*D(Jl)  )*TP(J5) 
19  CONTINUE 
DEL=0. 

no  15  1=1, MT 

J1=J  1+1 

CEL  = DEL  + D(  Jl  > 

15  CONTINUE 
J1=J1-MT 

SN=nEL 

DO  16  I=1,MT 
J5=J5+1 
J1=J  1+1 
S^=SN-D(  Jl) 

TP(  J5)=-D(  JD/OEL 
T( J5)=-( SN+,5#D( Jl ) ) *TP( J5) 
lb  CONTINUE 
Jl =J1-M1 
6 8 CCNTINUE 
3]2=  .5  + P! 

^r4  = NN 

DO  25  K=1,NPHI 
PH  = PI2>i‘(  X(K)  + 1.) 

PHT'=PH'‘FN 
SN=S  IN  ( .5=<'PH) 

CP(K  ) = 4.<'SN*SN 
R1  PI2*A(K ) 

cs  i i+cosc  phn; 

C 1 K )=  .5»C  P (K)*  ;c 
C2  K ) = C{iS(  PH)*C.> 
C3(K)=K1*SIN(PI )*S1N(PHN) 
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N?N=N2*'N 
N4N=  N2N=!=2 
OC  62  J=1,N4N 
Y( J) =0. 

Z ( J ) =0  . 

62  lCNTINUF 
U=(0.»1.) 
nr  59  j=i,NG 
FJ=FN/FS(J ) 

Ll=l 

L.’=2 

II  ( J .LF.MT  ) Ll=2 
IF(J.GT.NGM)  L2=l 
J9=( J-1 ) /MT 
JT-  J9-2) 

J5-  ( J9-2)*N2-1 
Sl=  1. 

DC  60  1=1, J 
I9=( I-l) /MT 
IT=I-frrT*( 19-2) 
J6=I9+J5 
Fl=FN/hS(l ) 
f;P  = F 5^(J)-f<S(I) 


ZP=ZS( J)-ZS(  I) 

f 2 = F.  p*p.p>2p#2p 

:p(I  .NF.J)  GC  TP  41 
i 1 ■=.  5 

F >.  0625*D(  J)-*'  M J) 

41  P:  =RS(  I )«RS  CJ) 

G1  = 0. 

G2=0. 

G3=0. 

G4=0. 

G5=0. 

G6  = 0. 

DP  6 1 K=1,NPHI 
R4=R2+F3*CF  (K) 

R5=SQrT(R4) 

Z1=S1/R5*(C(  S(R5)-U*SI  N(R5)  ) 

Y1=Z  l.  + U=  R5)/R4 

G1=C  1(K)*YH  G1 
G2=C2(K)*Yl-tG2 
G3=C3( K)#Y14G3 
G4  = r.  4(K)*2  14G4 
G5=C.2(K  )«Z1+G3 
G6=C  3(K)*Z1+G6 
61  CONTIFUE 

G3=U^'G3  „ 

Yl=(  RP»CV(J  )-ZP*S''(J))  *G2-RS(I)*CV(J)*G1 

Y2  = (RS(J)*  .V( I )*C  (J)-RS(I)*SV(J)«CV(I )-ZP*SV'  I)*SV( J) )*G3 
Y3=ZP«G3 

Y4=(  RP*CV<  I )-ZP*SV(I  ))  *G2+RS(  J I’t'CVI  I )*Gl 
Z1  = 'J’»(  SVt  1 ) ’I'SVf  J)*G5+C  V(  I ) *CV(  J)*G4) 

G1=-U*GA 
Z2=-SV( J)*G6 
G2=-FI*G4 
Z3=SV( 1 )*G6 


G3=F  3*04 


!i 
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Z4=U*(  liS-F  I*G3  ) 

Kl  = l 
K2  = 2 

IF(I.LE.MT)  Kl=2 
IF(I.GT.NGM)  K2=l 
DO  31  L=L1,L2 

lt=jt+mt*l 

J7=J6+L#N2 
DO  32  K=K1,K2 
KT=IT+MT*K 
TT=T(LH*T(KT) 

J1  = J7+K 
J2=J1+N 
J3=J1+N2N 
J4=J3+N 

Y( J1 )=TT*Y1+YJ  Jl) 

Y(J2  )=TT*Y2+Y(  J2) 

Y(J3 )=TT«Y3+Y( J3) 

Y(J4)=TT’f=Y4+Y(  J^) 

Z(  Jl  »=TT*Z1  + TP(LT)*TP{  .T)*G1  + Z(  JIJ 
Z( J2)=TT*Z2+TP( LT)*T(K  f)*G2  + 2(J2) 
Z(  J?)=TT*Z3+TP(KT)=<‘TILT)*G3^Z(  J3) 
Z(J4)=TT*Z4+Z(J4» 

32  CONTINUE 
31  CONTINUE 
6 0 CONTINUE 
59  CONTINUE 
N2P=N2+1 
KD1  = 1 
J1=0 
J5=0 

DO  11  J=1,N 
R1  = 0. 

DO  17  1=1, MT2 
J1=J  1+1 
J5=J  5+1 

Pl=F!l+T(  J1)*T(  J1)*PD(J5) 

17  CONTINUE 
J1=J1-MT2 
J5=J  5-MT2 
KD2=KD1+N 
KD3=KD1+N2N 
Kn4=KD3+N 
GUY  (KOI )-  i'  (KD4) 

Y(Kni)  =F1  + G1 
Y( KD2)=0. 

Y(KD3)=0. 

Y(KD4) =P1-G1 

Z(KD1)=Z(KD1)+Z(KDII 

Z(KD2)=Z(KD2»-Z(K03) 

Z(K03)=-Z(KD2) 

Z(K04)=Z(KU4»+Z(K04) 

IF(J-l)  26,27  26 
Zr  J1=J1+MT2 
j5sj  5+MT 
GO  TO  22 
26  KU1=KD1-1 
KU2»KD2-1 
KU3=Kn3-l 
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KU4=KD4-1 
KLl  = Kr)l-N2 
KL/?=Kn2-N2 
KL:1=KD3-N2 
KL*t=  KD4-N2 
R1=0. 

DO  18  1=1, MT 
Jl=  Jl+1 
J5=J5+1 

R1=R1+T(J1)#T(  J1-MT)*PD(  J5) 

I 8 continue 

J1=J  1+MT 
G1=Y(KU1)-Y(KL4) 

G2=Y (KU4)-Y(KL1) 

Y(KUl) =B1*G1 

Y(KU2) =Y(KU2)-Y(KL2) 

Y(KU3)=Y(KU3)-Y(KL3) 

Y(KUA)=R1+G2 

V (KLI)=R1-G2 

Y(KL2)=-Y(KC2) 

Y(KL3)=-Y(KU3) 

Y(KL^)=R1-G1 
Z(KU1)=Z(KU1»+Z(KL1) 
Z(KU2)=Z(KU2)-Z(Kl 3) 

Z(^U3)  =Z(KU3)-Z(Kl  2) 
Z(HJ4)=Z(KJA)+Z(KL4) 
Z(KL1)=Z(KU1) 

Z<KL2)=-Z{Kb3) 

Z(KL3)=-Z(KU2) 

Z(KL4) =Z(KU4) 

22  KD1=KD1+N2P 

II  continue 

IF(N.LT.3)  RETURN 
J2=N2 

DO  13  1=3, N 
J2=J2+N2 
Jl=I-2 
KL1  = I 

DP  lA  J=l,  J 1 
KU1=J2+J 
kU2  = KUl+N 
KU3=KJ1  + N2N 
KL4=KU3+N 
KL2  = KL  1 + N 
KL3=KL  1 + N2N 
KLA=KL3  + N 
Y(KL  1)=-Y(KUA) 

Y{KL2)=-Y(I'  U2  ) 

Y(KL3) =-Y(f  L3) 

Y(K.L4)=-Y(f  Ul) 

Z(KL1)=Z(KL  II 
Z(KL2)=-Z(KiJ3) 
Z(KL3)=-Z(KU2) 
Z(KL^)=Z{KU4) 

KL1=KL1  + N2 
l->  CONTINUE 
13  CP^’T1^UE 

return 

enp 
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III.  THE  SUBROUTINE  PLANE 


Description: 

The  subroutine  PLANE(NN,  N,  MT,  NT,  THR,  R)  puts  R**^  of  (49)  in 

nj 

R(j  + (m-l)*N+(L-l)*4*N)  where  j = l,2,...N,  nj=l  denotes  pq  = t6,  in=2 
denotes  pq  = i()9,  m=3  denotes  pq  = t<t>,  in=4  denotes  pq  = i|>4>,  ind  L*1,2,...NT. 
Here,  L denotes  the  Lth  value  of  < he  receiver  angle  0^.  The  first  5 
arguments  of  PLANE  are  input  variables.  Except  for  the  new  input  variable 
MT,  the  subroutine  PLANE  is  the  same  as  the  old  subroutine  PLANE  on  pages 
22-26  of  [10].  If  MT  = 2,  these  subroutines  are  exactly  the  same  as  far 
as  the  calculation  of  R in  terms  of  the  rest  of  the  input  variables  is 
concerned. 


The  integration  of  fj(t)  over  t inherent  in  (49)  is  approximate  i 
by  sampling  ^ j 2*MT  times  instead  of  4 times.  The  represen  ;ation  of 

pfj(t)  given  by  (66)  of  [9]  is  replaced  by  (77).  NN  is  the  value  of 
appearing  in  required  that  NN  > 0 but  this  requirement  ca  ises 

no  real  loss  of  generality  because  1?^  is  either  even  or  odd  in  n.  I is 
the  number  of  expansion  functions  lying  on  the  generating  curve.  Sp  :cifi- 
cally,  N Is  the  maximum  value  of  i in  (77).  THR(L)  is  the  Lth  vain  of 
tlie  receiver  angle  0^  where  L = 1,?,...NT.  The  variables  RS,  ZS,  S , CV, 
and  T appearing  in  the  common  statement  early  in  the  subroutine  PLA  iE  are 
input  variables  calculated  by  calling  the  subroutine  YZ  beforehand.  The 
calculated  values  of  these  variables  depend  only  on  the  second,  flf  h, 
sixth,  and  seventh  arguments  (NP,  MT,  RH,  ZH)  of  YZ. 

Minimum  allocations  are  given  by 


COMPLEX  R(4*NT*N) 

DIMENSION  7HR(NT),  BJ(M) 

COMMON  RS(NG),  ZS(NG),  SV(NG),  CV(NG),  T(2*MT*N) 

where 

NG  = (N+1)*MT 


and  M is  the  largest  of  the  /alues  of  M calculated  by  PLANE.  The  sug- 
gested allocation  BJ(50)  will  work  if  the  maximum  circumference  of  the 
body  of  revolution  is  less  tiian  26  wavelengths. 


«>8 


Most  of  the  statements  in  the  subroutine  PLANE  are  the  same  as 


or  very  similar  to  statements  in  the  old  subroutine  PLANE  listed  on 
pages  25-26  of  [10].  The  major  difference  is  in  the  calculation  of 
the  subscript  (IT+MT*K)  for  T and  the  subscript  J1  for  R.  Using 
reat  ming  similar  to  that  used  to  obtain  the  subscrip  :s  KT  and  J1  In  the 
sub  I Jutine  YZ,  we  arrive  at 


(IT  + MT*K)  = I + (19  + K-2)*M'’ 

J1  = 19  + K-1 

where 

19  = [('-!) /MT] 

The  above  Jl  is  valid  only  for  L=  1 . If  L 1,  then  (L-1)*A*N  must 
be  added  to  this  Jl. 


listing  nr  the  subroutine  plane 

SUBROUTINE  PLANE (NN»  N, MT,NT,ThR»P) 

COMPLEX  R(106^)  fUf  U1  ,U2 » R1  »K2 ,R3  ,114 
DIMENSION  THR(37j  , Bo  (50) 

COMf’ON  RSI  160)  f ZS(16  0)  ,SV(  1(  0)  fCV(l60)  tT(320 
NG=(  N+1)*^MT 
U=^  ( ( . , 1 . ) 

Ul  = 3 .141?  93»U**NN 

NA=4*N 

jr  = N4*NT 

Dl  2 2 J=1  JR 

R(  J) =0. 

22  Continue 
J5=-l 

DP  12  L=1,NT 
( s=r.c  (Thr (L) ) 

SN=2  SIN(  THR(L)  ) 

DO  13  ;=1»NG 
X=.25"^  S(I)*SN 
IF(X  .Lt..5P-7)  GO  TO  I 'i 
M=2. 8*X+13.-2./X 
IF(X.LT..5)  M=10.b+AL0ol0(X) 

IF(M.&E.(fvN  + 2)  ) GO  TO  19 

18  BJ1  = 0. 

BJ2=0. 

BJ3=0. 

IF(NN.Ft.l)  bJl=l. 

IECNN.EQ.O)  8J2*1. 

GO  TO  24 

19  Bj(M)  = ' . 

JM=M-1 

BJ(JM)=1. 
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DO  14  J*3fM 

BJ(  JM)  =JM/X*BJIJM»1)-BJ(  JM+2) 

14  CrNTlNUE 
S*0. 

DO  I 5 J = 3,M,2 
S=S*BJ(  J) 

I 5 CONTINUE 

S=BJ (1 l♦2.*S 
BJ2*BJ  <NN>1)  /S 
BJ3=  BJ  (NN+2  ) /S 
J3 

IFINN.GT.O)  BJ1=BJ(NN)/S 
24  APG=ZS(I)*CS 

U2  = U1*(CCS«  ARG)  + U»SIN(  ARG) ) 
k4=( BJ3-BJ1 )*U*U2 
R2  = ( BJ3+Bjn*U2 

R 1=-BJ2*CV(  I )*SN4<U2+CS*SV(  I ) *R4 

R3*SV( I )*R2 

E2=-CS>f'R2 

I9=(  I-l  )/MT 

IT=I+MT*(I9-2) 

J^=I9+J5 


IS.  -c. 

11(19. EQ.O)  Kl*2 
IE!I9.EQ.N)  K2*l 
OL  20  K = K1,K2 
TT=T(IT+yT*K) 
Jl*J7+K 
J2*Jl+N 
J3*J  2+N 
J4*J3+N 

R(J1 )=TT*R1+R( Jl) 
R(  J2  1=TT*R2+R(J2) 
R( J3 )=TT*P3+R( J3) 
R( J4 )=TT*R4+R( J4) 

20  continue 

13  C ^NT  INUE 
J 5*J  5+  '(4 

12  continue 

RETURN 

END 
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IV.  THE  SUBROUTINES  DECOMP  AND  S< LVE 


Description: 

The  subroutines  DECOMP (N, I PS, UL)  and  SOLVE (N, IPS, UL.B.X)  solve 
a system  of  N linear  equations  in  N unknowns.  These  subroutines  will 
be  used  in  Section  V to  solve  the  matrix  equation  (39).  The  input  to 
DECOM*  consists  of  N and  the  N by  N matrix  of  coefficients  on  the  left- 
hand  side  of  the  matrix  equation  stored  by  columns  in  UL.  The  output 
from  DECOMP  is  IPS  and  UL.  This  output  is  fed  into  SOLVE.  TliU  rest  of 
the  input  to  SOLVE  consists  of  N and  the  column  of  coefficients  on  the 
right-hand  side  of  the  matrix  equation  stored  in  B.  SOLVE  puts  the 
solution  '.o  the  matrix  equation  in  X. 

Minimum  allocations  are  gl\'en  by 

COMPLEX  UL(N*N) 

DIMENSION  SCL(N),  IP5(N) 

in  DE  ■;OMP  and  by 

COMPLEX  UL(N*N),  B(N),  X(N) 

DIMENSION  IPS(N) 

in  SOLVE. 

More  detail  concerning  DECOMP  and  SOLVE  is  on  pages  46-49  of  [13] 

LISTINd  CF  THE  SUBROUTINES  DECOMP  AND  SOLVE 

SLBpnUTiNE  DFCOMPt Nt IPS.UL) 

COMPLEX  UL ( 3136) fPIVOT.EM 
DIMENSION  SCL(  56),IPS(5t  ) 

DC  5 I=1.N 
IPS(  I)=I 
RN  = 0. 

JUI 

Dr  ? j = i,N 

ULM=A6S(REAL<UL( J1 ) ) )+ABS(AI MAG(UL (Jl) ) ) 

J1=J 1+N 

IF(RN-ULM)  If2f2 

1 RN=UL  1 

2 CCNTi;  UE 
SCL(I  =1./RN 

5 CfiMIi  UF 
NM1  = N-1 
K2  = 0 

00  17  K=1,NM1 
BIG=0. 

DO  11  I*KtN 
IP=^I  PS  ( I ) 

1PK=IP+K2 
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f 
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S1ZE=( APS( RFAL(UL( IPK)  ) ) ^ABS ( Al MAG ( UL ( IPK) ) ) )*SCL( I'  ) 
IF(SIZE-BIG)  lltllflO 
K BIG=SIZF 
IPV=  I 

11  CONTINUE 

IFUPV-K)  l<^,15tlA 
14  J=IPS(  K) 

IPS<K)=IPS( IPV) 

IPS( IPV)=J 
1 5 KPP=IPS(  K)+K2 
PIV0T=UL(KPP) 

KP1  = K+1 

DO  16  I =KP1 tN 

KP=KPP 

IP=IPS( I )+K2 
EM=-UH  IP) /PIVOT 
18  UL(I  P)=-EM 
DO  16  J=KP1  ,N 
lP=rIP+N 
KP=K  P+N 

UL(  P)=UL(IP)  + EM*ULIKP) 

1..  cor''  INUE 
K2=K?+N 
1 CONTINUE 
RETURN 
END 

SULPOUTINE  SOLVEINtI PStULtBfX) 

COMPLEX  UL (3136) ,B( 56) ,X( 56) ,SUM 
DIMENSION  IPS(5b) 

NP1=N+1 
IP  = IPS(  1 ) 

<(1)  =B(  IP) 

I'O  2 I=2tN 
IP  = IPS(  I) 
ipe=  IP 
JM1=I-1 
SUM=0. 

DC  1 J*1,IM1 
SUM=SUM+UL( IPI^XIJ) 

1 IP=IP'^N 

2 X( I )=P( IPB)-SUM 
K2  = N*(N-1) 

IP=I PS (N)+K2 
X(N) =X(N)/UL( I P) 

Df  4 IBACK=2,N 
I=NP1-I PACK 
K2=k2-N 
rpi=iPS(  n+K2 
.PI*  i+i 
: UM*o. 
ip=i  PI 

DO  3 J=IP1 »N 
IP=I P+N 

3 SUM*SUM+UL ( IP)*X(J) 

4 X(I)=|X(I)-SUM)/ULIIPI) 

RETURN 

END 
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V.  THE  MAIN  PROGRAM 

Description: 

The  main  program  calculates  the  electric  and  m.  gnetic  currt  ts 
(58)-(59)  and  the  normalized  scattering  cross  sections-^  and  ' /O) 

A A 

for  the  6 polarized  axially  incident  (6  = 180°)  plane  wave  (50)-(51). 

The  main  program  calls  the  subroutines  YZ,  PLANE,  DECOMP,  and  SOLVE. 

The  main  program  is  short  and  simple.  It  is  a representative  applica- 
tion of  the  theory  in  Part  One  of  this  report. 

Input  data  is  read  early  in  the  main  program  according  to 

READ(1,10)  NP,  NPHI,  MT,  NT 

10  FORMAT (41 3)  ’ 

READ (1,11)  BK,  UR,  E’l,  ALP,  BET 

11  FORMAT(5E14.7) 

READ(1,12) (RH(1),  I=L,  NP) 

READ(1,12)(ZH(I),  1=  , NP) 

12  FORMAT(10F8.4) 

READ(l,ll)(xrK),  K=l,  NPHI) 

READ(1,11)(A.K),  K=l,  NPHI) 

The  input  variables  NP,  NPHI,  MT,  RH,  ZH,  X,  and  A are  very 
similar  to  variables  of  the  same  names  in  the  argument  lis»  of  the  sub- 
routine YZ.  In  summary,  (RH(i),  ZH(i)),  i=l,2,...NP,  are  the  c lindrlcal  1 

coordinates  (Pj^»  z^)  on  the  generating  curve,  2*MT  is  the  numbe  of  values 

of  t at  which  f.(t)  is  sampled  for  the  purpose  of  numerical  integration,  i 

^ 1 
and  X and  A are  respectively  the  abscissas  and  weights  for  the  NPHI  point 

Gaussian  luadrature  Integration  in  iji.  i 

Tht  scattering  cross  sections  are  evaluated  at  receiver  ai  gles 

0^  = (J-1) I / (NT-1)  radians  for  J=1,2,...NT.  BK  is  the  propagation 

constant  k in  the  external  medium.  This  k appears  in  \.50)-(51).  UR 

bd 

and  ER  are  rjspectively  the  relative  perme  ibility  — and  the  relative 

Ed  *^6 

permittivity  — of  the  body  of  revolution.  Here,  p , and  e , are  the 
c d d 

e 

permeability  and  permittivity  of  the  (diffracting)  body  of  revolution  | 
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and  p and  e are  those  of  the  external  medium.  ALI  and  BET  are  r^  spec- 
e e 

tively  thi!  constants  a and  6 appearing  in  (39).  Tht  PMCHW  solutlo  is 
obtained  by  setting  a = 6 = 1.  The  Muller  solution  is  obt<  ined  if  a and 
6 are  given  by  (22)  and  (23). 

Minimum  allocations  are  given  by 

COMPLEX  YE(4*N-N),  ZE(4*N*N),  R(4*NT*N) , 

B(4*N),  YD(4*N*N),  ZD(4*N*N) , Y(16*N*N), 

C(4*N) 

DIMENSION  RH(NP),  ZH(NP),  X(NPHI),  A(NPHI), 

THR(NT),  RC(N),  IPS(4*N) 


where 


N = [(NP-1)/MT]  - 1 

Statement  38  puts  Y^^  defined  in  accordance  with  (75)  oy 


le 


tt 

le 

le 

yH 

le 

le 

(84) 


in  YE  and  defined  in  accordance  with  (76)  by 


"le 


le 


At 


"le 


.A  ■ 


(8') 


in  ZE, 


"le  lej 

Storage  of  Y^^^  and  Z^^  is  by  columns. 
Stateiitent  39  puts  the  matrix  Rj^  defined  by 


ite 


r1 


(86  1 


[ I'l  -1 

in  R.  The  column  vectors  on  the  right-hand  side  of  (86)  appear  in  (56). 
For  receiver  angle  0^  = THR( F) , R^  of  (86)  is  stored  by  columns  in 
R((J-1)*4*N+1)  through  R((J*<*N). 
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DO  loop  22  uses  (56)  to  store  the  right-hand  side  of  (39)  i B. 

DO  loop  22  also  puts  2/ in  RC  where  is  the  cylii  Irical 

coordinate  radius  evaluated  at  the  peak  of  the  triangle  function  ■ nherent 
in  fj.(t)  . 


Statement  40  puts  defined  by 


Id 


in  YD  and  Z, , defined  by 
id 


^Id 


tl 

Yt-O 

l.l 

Id 

(fit 

Id 

’^Id 

tt 

Id 

^Id 

Id 

^Id 

in  ZD.  Storage  of  Y^^  and  is  by  columns. 


(87) 


(88) 


Nested  DO  loops  26  and  27  >ut  the  first  two  columns  of  submatrices 
on  the  left-hand  side  of  (39)  in  Y.  The  index  J of  DO  loop  26  denotes 
the  Jth  column  of  the  composite  square  matrix  on  the  left-hand  side  of 
(39). 


Nested  DO  loops  23  and  29  put  the  third  and  fourth  columns  of  sub- 
matrices on  I he  left-hand  side  of  (39)  in  Y.  The  index  J of  DO  loop  28 
denotes  the  (2*N+.J)th  column  of  the  composite  square  matrix  on  the  left- 
hand  side  of  (39). 

Statements  41  and  42  solve  the  matrix  equation  (39)  for  the  com- 
posite column  vector  consisting  of  V^,  V^,  and  I^.  This  composite 
column  vector  is  stored  in  C. 


At  the  peak  of  the  Jth  triangle  function,  the  n=J  term  of  the 
equivalent  electric  current  (58)  reduces  to 
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DO  loop  31  prints  th(?  real  and  imaginary  parts  of 


U1  j , 

the  real  and  imaginary  parts  of 

and  the  magnitudes  of  U1  and  U2.  Here,  U1  is  the  t component  of  the 


equivalent  electric  currt  nt  in  the  E plane  and  U2  is  the  ij)  component 
of  the  equivalent  electric  current  in  the  H plane  when  the  y componei  t 
of  the  incident  magnetic  field  is  minus  one  at  the  origin. 

At  the  peak  of  the  Jth  triangle  function,  the  n= : term  of  thi 
equivalent  magnetic  current  (59)  reduces  to 

DO  loop  34  prints  the  real  and  imaginary  parts  of 

B(J)  - (2J/(k.-T.j^p)v“ 
and  the  real  and  imaginary  parts  of 

and  the  magnitudes  of  B(J)  and  U2.  Here,  B(J)  is  the  t co  iponent  of 
the  equivalent  magnetic  current  in  the  H plane  and  U2  is  the  iji  component 
of  the  equivalent  magnetic  current  in  the  E plane  when  the  x component 
of  the  inc  dent  electric  field  is  minus  one  at  the  origin. 

Nested  DO  loops  35  aid  36  calculate  and  print  and  (70). 

Inner  DO  loop  36  accumulates  the  portion 
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-wm 


^ ^ -te-te  ^ 


of  (62)  in  ET  and  the  portion 


-tentte  _ ^<t)9^(|)e  st(j)>te 


1 '1  ^1  ^1 


of  (63)  in  EP.  The  W1  ar J W2  printed  in  DO  loop  35  are  respectively 
0 0 

and  -#  for  receiver  ingle  6 = ( J-1) *7r/ (NT-1)  radians. 

A A IT 

Suggested  modifications  of  the  main  program  are: 

1)  Changing  Che  normalization  of  the  scattering  patterns. 

For  example,  one  could  replace  (70)  by  (69).  All  scat- 
tering patterns  in  Part  One,  Section  V are  plots  of  (69). 

2)  Removing  the  restriction  th-.  t the  values  of  the  input 

arguments  NP,  NPHI , , RH,  ZH,  X,  and  A of  the  subroutine 

YZ  be  the  same  in  call  statements  38  and  40.  This  modifi- 
cation is  indicated  by  (74). 

3)  Generalizing  from  axial  plane  wave  incidence  to  oblique 
plane  wave  incidence. 

The  above  three  modifications  can  be  realized  without  tampering  with 
any  of  the  subroutines  YZ,  PLANE,  DECOMP,  and  SOLVE. 

The  sample  input  and  output  accompanying  the  listing  of  the 


main  program  is  for  the  dielectric  sphere  with  ka  = 1 and  f 
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LISTING  f'F  THF  MAIN  PROGRAM 


//PSM  JOB  (XXXX,XXXX,1»2) ,*MAUT2,J  IE*  »RFGI0N=200K 
//  EXEC  tvATFIV 
//GD.SYSIN  00  * 

*JOP  MAUTZ,TIME=lt  PAGES=40 

C SUBROUTINES  YZ,  PLANE,  DECOMP,  AND  SOLVE  ARE  CALLED. 

COMPLEX  YE(78A) ,ZE(784),R(106A) ,B(56) ,Y0(784),Zn(784l  Y(313  > 
CCMPLFX  C( 56) ,U,U1,U2,ET,EP,C0NJG 

01  ME  NS  ION  RH(lbl),ZH(161),X(48),A(48)  ,THR(37  ) ,RC  ( 19)  .IPS(5<;) 
READ(1,10)  NP,NPH1 ,MT,NT 

10  FURMAT(4I3) 

READ(1,11)  BK,LK,ER,AL  P,BET 

11  FORM  ATI 5E14.7) 

READ(1,12)(RH(I),I^1,NP) 

READd  ,12)  ( ZH(  I)  ,1  = 1 ,NP) 

12  FORMAT!  10F8. 4) 

READd  , 1 1)  ( X(  K)  ,K=1,  NPHl  ) 

RFADd,  11)(A<K),K=1,NPHI) 

WRITEI3,13)  NP,NPHI,MT,NT 

13  FORMAT!'  NP  NPHI  MT  NT*  / ! 1 X . 13 , 15 , 2 13  ) ) 

WF1TF!3,14)  BK,UR,ER,ALP,BET 

14  FORMAT!  7 Xt'BK*  , 1 2X  , * UR  * , 12X  , ' ER  * , 1 1 X,  • AL  P'  , 1 1 X,  • B ET*  / ! 1 X ,5  F 4.7)  ) 
WRITF13,15)!RH(I),I=1,NP) 

15  FOPr  AT(«  PH' /(  1X,10F8.4)  ) 
kvRnF!3,16)  !ZH!I),I=1,NP) 

16  FORMAT!'  ZH' /<  1X,10F8.4)  ) 
k^RIT  F!3 ,17)  (X(  K)  ,K=1  ,NPHI  ) 

17  FORMAT!'  X' / ! 1 X , 5E14 . 7 ) ) 

WRITE!  3,18)  ! AIK)  ,K=1,NPHI) 

18  FORMAT!'  A ' / ! IX , 5E  14 .7  ) ) 

DO  19  J=1,NP 

RH( J )=BK=*RH(J  ) 

Zh! J)=BK*ZH!J) 

19  continue 

38  CALI  YZ(1,NP,NPH1 ,-l,MT,RH,ZH,X,A,YL,ZE) 

WRnF!3,20)  !YF!I),I=l,4),(ZF!i)  ,1=1,4) 

20  FORMAT!'  SOME  ELEMENTS  OF  YE  AND  ZE ' / ! IX , 4E 14. 7) ) 

PI =3 .141593 

nT  = P I/INT-1 ) 

PC  21  J=1,NT 
THF  ! J)  =!  J-1 1’^DT 

21  CCNTINUF 
N=!NP-l) /MT-1 

39  CALL  PLANE! 1,N,MT, NT, THR,R) 

N2=2«N 

N3=3*N 

N4=4=N 

NTr.  = !NT-l)«N4 
no  22  l=ltN 
J3=I+NTN 
E!I)=R!  J3) 

B!  I'-N)=-R!  J3+M  ) 

B!  I + N2  ) = -R(J3  + N:-  ) 

B!  1 + N3)=-R!  J3  + N2  ) 

RC  (I  ) = 2./RH(MT*d  1) 

22  Ct  JT  INUE 

WRITE!  3,23)d  ( I)  ,1=1  ,N2) 

2?  FOFMAT!'  HALF  OF  'HE  ELEMENTS  OF  B' / ( 1 X ,4E  14 . 7 ) ) 

FMf  ,CRT!  UR*ER) 
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DO  24  I=1»NP 
RH(  I )=EM’!‘RH(  I ) 

ZH(I ) = EM*2H(  I) 

24  CDImTINUF 

40  CALL  YZ(ltNP,NFHIflfMT,RHtZH,X,A,''r),ZD) 

WRITE(  3,25)  (YD(I),I=1,4)  ,(ZD(I)  *1=^  lt4) 

25  FCRMATC  SOME  ELEMENTS  OF  YO  AnD  D • / ( IX , 4E14. 7)  ) 
n*SQRT(UR/EP) 

ALPD=ALP*D 
faETD=BET/D 
JY=0 
J1  = 0 

DC’  2 6 J-1,N2 
DP  27  1=1, N 
JY=JY+1 
J1=J1+1 
J2=J1+N 

Y(JY )=YE( J2 )+ALP*YD( J2) 

. Y(  JY+N)=-YF  ( J1  )-ALP»YD<Jl) 

Y( JY+N2)=ZE(J2)+BET0*ZD( J2) 

Y(JY+N3)=-ZE( J1)-DETD* ZOCjl  ) 

2 7 CCNTINUF 
JY=JY+N3 
J1=J  1+N 

26  CGNTINUF 
J1=0 

DO  28  J=1,N2 
DC  29  1=1, N2 
JY=JY+  1 
J1=J  1+1 

Y(  JY)  = ZE(J1  )+ALPD’»ZD(J  1) 

Y( JY+N2)=YF ( J1 )+BET*YO( Jl) 

29  cr  TINUE 
JY^JY+N2 

2 8 CONTINUF 

41  CALL  DECCHP  (N't  ,I  PS,Yi 

42  CALL  SLLVE(N4,IPS,Y, j,C) 

V<R  1TE(  3,30) 

3 0 FCRMATCO  PEAL  JT  IMAG  JT  REAL  JP  IMAG  JP 

IMAG  JP*  ) 

U=(0.,1.) 

DC  31  J = 1,N 
RP=RC( J> 

ET=U*kR 
B( J)=ET+C( J) 

J1=J+N 

B(  J1  ) = RP.«C(  Jl» 

U1=RP*C ( J+N2) 

U2=ET*C( J+N3) 

W1=CABS(  LI  ) 

K2=CABS(U2) 

k»PITE(3,32)  U1,U2,W1,W2 

32  FORMAK  1 X,6E11.4) 

3 1 CONTINUf 

)»FITE(3.33) 

33  FfrMAT(‘0  REAL  MT  IMAG  IT  REAL  MP  IMAG  MP 
IMAL  MP*J 

no  34  J=1,N 
U2=3(J+N| 
wi  = :ABS(B(jn 


W?=CABS(U2I 

ksPITE(3f32)  B(  J)  ,U2,k«l  ,W2 
3 4 c.ONTINUF 
WPIT  F(3,37) 

37  FCRMAK’O  SIGTHETA  SIGPHI') 
CCiN=  .2  5 /F  1**3 
JF1  = 0 

OP  35  J=1,NT 
FT=0. 

FP=0  . 

DP  36  I=ltN2 
JC2= I+N2 


Jf  1 = JP  1 + 1 
JP2=  JP  1*P'2 

rT=P7+K ( JR2 )*C(I)+R4  JR1)*C( JC2) 
FP=FP-P,(  JRl  )*C(  I l+R(  JR2)*C(JC2) 

36  CPMINuE 
JP.1=JP.1+N2 

CN*FT*tONJG(  ET ) 

W2=f  CN*FP*CCNJG( EP ) 

ViRITF(3,32)  Wl,W2 
3 5 (PNTinuE 
STOP 
FOP 

$DAT  A 
21  20  2 
C.l OPOOOOF+Ol  0 

0.  00  00  0.  1564 

1.  0000  C. 51877 

0.  0000 

-1.  0000  -0.9877  -0.9511  -0.  8910  -0.  8090  -0 

0. 0000  0.  1564  C.3090  0.4540  0.5878  0 

1.  0000 

-0.9931 286 F+  00-0.9639719E+00-0.9122344F+00- 
-0.6  3605‘’7F  + 00-0. 510  8670F+  00-0.373  7061E  + 0 0- 
0.7652692  F-01  0 . 2 2778  59E+ 00  0.  373706  1 E+00 
0.7463316F+01)  0. 83911  70E^C0  0.912  2344  8+00 
0.1 76 1401E-01  0.40601438-01  0.6267205E-01 
0.1U1945F+00  0. 1316886E+00  0. 1420961 E+00 
0.1527534F+00  0 . 149 1 730E+ 00  0. 142 096 1 F+0 0 
0.1 019301  E + 00  0. 8327674E-01  0 .626 720 5 F-0 1 
$ST3P 
/* 

// 


lCOOOOOF+01  0.4000000F+01 
0.3090  0.4540  0.5878  0 

0.9511  0.  8910  0.8090  0 


0. I OOOOOOE+01 
. 7071  0.8090 
.7071  C.587P 

.7071  -0.5878 
.7071  0.809C 

9.8391170E+0C- 
0.2277859E+0r- 
0.  5108670E+C0 
0.9639719E+00 
0. 8327674E-01 
0. 1491730E+00 
0. 1316886E+00 
0.4060143E-01 


O.lOOOOCOF+01 
0.8910  0.9511 

0.4540  0.3090 


-0.4540 

0,8910 


0.3090 

0.0511 


■0.746331  )E+00 
0.76526528-01 
0.63605  37 E+OO 
0.9931286F+00 
0.10193C1E+00 
0.1527534E+00 
0.11819^58+00 
0.17614i;lE-01 


0.9877 
0. 1564 

-0. 1564 
0.9877 


PP  INTF  0 ( UTPUT 
NP  NPMl  MT  NT 
21  20  2 IS 

8K  UP  FP  ALP  BET 

0.1  OOlOf  0F+ 91  0.  lOCCOOOL  + 01  0.400  0000^  + 91  0 . lUOOOOOE+01  0 . 1 9000 00 F+ 1 
PH 


j 0.  00  00 

0.1 564 

0.3090 
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0.  7071 
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1.  00  00 

0.9877 

0.9511 

0.  8910 

0. 8090 

0. 7071 

0.5878 

0.4540 

0.3090 

0.  1564 

0.  ocoo 
ZH 

j - 1 . or  00 

-0.  98  77 

-0.95 11 

-0.  8910 

-0.8090 

-0.  7071 

-0.  5878 

-0.4540 

-0.3090 

-0. 1564 

C.  0000 

C. 1564 

0.3090 

0.  4 .40 

0.59  78 

0.  7071 

0.  8090 

0.8910 

0.9511 

0.9877 

1.  0000 
X 
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0.49421258-01-0. 13746378+00 
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0.14658608+00  0.19117348+00 


FFAL  I 

■0.63988  + 0 0—  0 . 
-0.8349F+00-0. 
-0.1C258+0  1-0. 
-0.11278  + 01-0. 
■0.1094F  + 01  0. 
•0. 87318  + 00  0. 
-0.541 OF+00  0. 
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RE  AL  M T I 
-0.1  Of  4F  + 00-0. 
-0.22238  + 00-0. 
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-0.47318+00-0. 
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REAL  MP 
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RO.IECrr  SUMMARY  AND  RliCOMMEND  TIONS 


I.  SUMMARY  OF  WORK  PERFORMED 

The  purpose  of  this  project  was  to  inves  igate  applic  tions  of 
the  method  of  moments  to  obtain  vector  solutions  of  integro-  lifferentia] 
equations  Df  radiation  and  scattering  problems  for  se  ^eral  types  of 
boundary  conditions.  This  first  required  a study  of  techniques  which 
apply  to  the  general  formulation  of  matrix  solutions  to  multiple  regions 
ard  multiple  sources.  After  the  gent^ral  f(  emulation  was  complete,  the 
S(  lution  was  to  be  applied  to  repres  ntative  problems  such  as  arrays  of 
cavity-backed  and  waveguide-fed  aper  ures. 

Thu  first  problem  was  to  appl  the  techniques  to  aperture  antennas 
in  conducting  bodies  of  revolution.  This  involved  the  use  of  previously 
developed  programs  for  radiation  and  scattering  from  conducting  bodies  of 
revolution  [14,15].  These  iirograms  are  a solution  to  the  E-field  integral 
equation,  and  it  is  known  tliat  solutions  to  this  equation  are  not  unique 
at  certain  cavity  resonant  frequencies  [9].  The  alternative  approach  of 
using  the  H-field  integral  equation  also  has  solutions  that  are  not  unic  le 
at  the  same  cavity  resonant  frequencies  [9].  One  method  to  obtain  an 
Integral  equation  which  does  have  unique  solutions  at  the  cavity  resona  t 
frequencies  is  to  take  a linear  combination  of  the  E-field  and  H-field 
equations  [1].  Such  a procedure  leads  to  the  combined-field  integral 
equation.  However,  no  computer  programs  were  available  for  solutions  to 
the  combined-field  equation  for  bodies  of  revolution,  so  the  first  task 
was  to  develop  such  programs.  This  resulted  in  our  firs  two  Interim 
Technical  Reports  [9,  10].  Abstracts  of  these  reports  aie  given  ir 
Section  II  below. 

[14]  J.  R.  Mautz  and  R.  F.  Harrington,  "Generalized  Network  Parameters 
for  Bodies  of  Revolution,"  Report  AFCRL-68-0282 , Contract  No. 
F19628-67-C-0233,  Air  Force  Cambridge  Research  Laboratories,  May  1968. 

[15]  R.  F.  Hai rington  and  J.  R.  Mautz,  "Radiation  and  Scattering  from  Bodies 
of  Revolution,"  Report  AFC  . L -69-0305,  Contract  No.  F19628-67-C-0233, 

Air  Force  Cambridge  Research  Laboratories,  July  1969. 
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Development  of  the  combii  ed-field  sol  ition  for  conducting  bodies  ■ f 
revolution  required  the  evaluati  m of  both  t le  E-fleld  and  the  H-field 
integral  operators.  These  two  operators  are  also  those  required  to  solv  ■ 
the  problem  of  electromagnetic  scattering  by  dielectric  bodies.  Hence, 
our  next  task  was  to  incorporate  the  computer  programs  for  E-field  and 
H-field  operators  into  programs  for  dielectric  bodies.  The  theory  of  th  s 
work  constitutes  Part  One  of  this  report.  The  computer  programs  and  dot  i- 
mentation  of  them  are  given  in  Part  Two  of  this  report. 

The  next  step  in  the  project  was  to  a )ply  these  st'lutions  to  the 
problem  of  waveguide-fed  and  cavity-backed  apertures  in  onducting  bodies. 

These  waveguides  and  cavities  may  or  may  not  be  dielectr  c filled.  Pre- 
liminary work  for  waveguide-fed  apertures  was  given  in  Ri  port  No.  12  of 
the  preceding  project  [16].  Application  of  the  theory  to  arrays  of  waveguide- 
fed  and  cavity-backed  antennas  li  conducting  planes  is  given  in  another  two 
reports  [17,18].  It  was  planned  to  use  a similar  type  of  solution  for  aj er- 
tures  in  a conducting  body  of  re/olutlon.  This  solution  would  use  the  computer 
programs  given  in  the  present  report  (Part  Two)  to  determine  ;he  apertur  ■ 
admittance  matrix.  However,  further  worl  must  be  done  before  this  solut  ion 
can  be  made  available. 

1 1 . INTERIM  ECHNICAL  REPORTS 

during  the  project  two  Interim  Technical  Reports  were  written  and 
published  by  RADC.  These  reports  and  their  abstracts  are  listed  below. 

1.  "H-field,  E-fi  'Id,  and  Combined -fie Id  Solutions  for  Bodies  of 
Revolution,"  by  J.  I . Maut z and  R.  F.  Harrington,  Interim  Technical  Report 
RADC-TR-77-109,  Marc  i 1977. 

[16]  I.  R.  ffautz  and  R.  F.  Harrington,  "Transmission  fron  a Rectangular 
Wavegui Je  into  Half  Space  through  a Rectangular  Aperture,"  Report 
RADC-TR-76-264,  Rome  Air  Development  Center,  August  1976. 

[17]  J.  Luzwick  and  R.  F.  Harrington,  "A  Peactively  Loaded  Aperture  Antenna 
A’ ray,"  Report  TR-76-10,  C< ntract  No.  N00014-76-C-0225 , Office  of  Naval 
R' search,  September  1976. 

[18]  J Luzwick  and  R.  F.  Harrington,  "A  Solution  for  a Wide  Aperture 

R»  actively  Loaded  Antenna  Array,"  Report  TR-77-1,  Contract  No.  N00014- 
7(-C-0025,  Office  of  Naval  Research,  January  1977. 
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Abstract!  H-field,  E -field,  and  combined-field  solutions  are 
developed  for  the  electric  s irface  currei  t and  far  scattered  fields  of 
a perfectly  conducting  body  of  revolutioi  excited  by  an  incident  plane 
wave.  These  solutions  are  obtained  by  applying  the  method  of  moments 
to  the  H-field,  E-field,  and  combired -field  integral  equations,  respec- 
tively. The  H-field  integral  equation  is  obtained  by  requiring  the 
tangential  magnetic  field  to  be  zero  just  inside  the  surface  S of  the 
body  of  revolution.  The  E-fleld  integral  equat ion  is  obtained  by  re- 
quiring the  tangential  electric  field  to  be  ze  o on  S.  The  combined 
field  integral  equation  is  a linear  combination  of  the  H-fleld  and  E- 
field  integral  equations.  Computations  show  that  both  the  H-field  an» 
the  E-field  solutions  deteriorate  near  Internal  resonances  of  the  con- 
ducting surface  S,  but  that  the  combired  field  solution  does  not. 

2.  "Computer  Program.';  for  H-fjeld,  I -field,  and  Combined-fie’ d 
Solutions  or  Bodies  of  Revolution,"  by  J.  K.  Mautz  and  R.  F.  Harrir  ;ton, 
Intt  rim  Te  hnical  Report,  RADC-TR-77-215 , June  1977. 

Abst  ract:  A computer  program  is  given  to  implement  the  H-fl  Id, 
E-field,  and  Combined-field  solutions  given  in  Interim  Technical  Report 
Ri\l)  :-TR-77-109  for  a perfectly  condiu  ting  body  of  revolution  excited  by 
an  )blique  plane  wave  Incident  field  The  program  consists  of  s-everal 
sub  outines  and  a main  program.  ITie  main  program  calcula.es  thi  elec- 
trit  current  on  the  b')dy  of  revolution  and  the  blstatic  scattering  cross 
section  per  square  wavelength.  Some  examples  of  computations  are  tiven 
and  discussed. 

III.  JOURIAL  PUBLICATIONS 

During  the  project  the  following  papers  on  work  related  to  tiie 
project  have  been  publlsiied  or  accepted  by  technical  journals. 

1.  R.  F.  Harrington  and  J.  R.  Mautz,  "A  Generalized  Network 
Fornulatlon  for  Aperture  Problems,"  IEEE  Trans . . vol.  AP-24,  No  6, 
Noveirber  1976,  pp.  870-873. 

2.  R.  F.  Harrington  and  I.  R.  Mnutz,  "Electromagnetic 
Transmission  through  an  Aperture  In  a Conducting  Plane,"  AEU,  vol.  i. 

No.  2,  February  1977,  pp.  81-87. 
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3.  R.  F.  Harrington  and  J.  R.  Mautz,  "Conqjutatioi  il  Methods  for 
Transmission  of  Waves  through  Apertures,"  to  appear  in  E ectromagnetic 
1 Scattering,  edited  by  P.L.E.  Uslenghi,  Academic  Press,  1977. 

1 A.  Yu  Chang  and  R.  1'.  Harrington,  " v Surface  Formulation  for 

( 

!'  Characteristic  Modes  of  Material  liodies,"  accepted  for  publicatl  m in 

! the  IEEE  Transactions  on  Antennas  and  Propagation. 

I 5.  J.  L.  Luzwick  and  R.  F.  Harrington,  "A  Reactively  Loaded 

[ Aperture  Antenna  Array,"  accepted  for  publication  in  the  IEEE 

i Transactions  on  Antennas  and  Propagation. 

i 

^ I\ . RECOMMENDATIONS  FOR  FUTURE  WORK 

On  the  basis  of  results  obtained  during  the  present  project  'le 
following  topics  appear  promising  for  future  research. 

1.  A computer  program  should  be  developed  for  computing  the 
electromagnetic  behavior  of  waveguide-fed  and  cavity-backed  apet cures 
in  bodies  of  revolution.  The  waveguides  could  be  either  air  filled 
or  dielectric  filled.  The  program  would  make  use  of  the  matrices 
evaluated  for  dielectric  bodies  of  revolution,  appearing  in  this  report. 

2.  A combined-source  !■  olution  should  be  develope  i for  conducting 
bodies  of  revolution,  similar  o that  used  by  Bolomey  ai  1 Tabbara  [19] 
for  two-dimensional  scattering  problems.  The  advantage  of  the  combined- 
source  forinul.ition  is  that  it  a ^sures  uniqueness  of  the  solution  for  the 
E-field  (or  H-field)  at  all  frequencies.  It  is  probably  related  to  the 
combined-field  formulation,  given  in  the  first  Interim  Technical  Report. 

3.  Combined- f ield  and  combined-source  solutions  should  be  in- 
vestigated for  aperture  problems.  With  the  formulation  presently  used, 
the  admittance  matrix  for  a cavity  problem  becomes  singular  at  resonant 
frequencies  of  the  cavity.  It  is  Celt  tiiat  a combined-field  or  co  ibined- 
souroe  formulation  would  overcome  this  difficulty. 

[19]  .1.  C.  Bolomey  and  W.  la^bara,  "Numerical  Aspects  on  Coupling  between 
Complimentary  Boundary  .alue  Problems,"  IEEE  Trans.,  vol.  AP-21, 

No.  3,  May  1973,  pp.  356-363. 
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4.  General  programs  of  the  type  deve]oped  on  this  project  for 
bodies  of  revolution  should  be  developed  for  other  kinds  of  bodies. 
While  it  is  probably  too  dil ficult  at  this  time  to  develop  programs 
for  arbitrarily  shaped  bodies,  it  is  possible  to  develop  them  for  oth>r 
classes  of  shapes,  such  as  cylinders  of  arbitrary  cross  section,  elthf  r 
infinite  or  finite  in  length. 
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